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I. INTRODUCTION 



\I-ary differential phase shift keyed (M-DPSK) modulation is a widely used 
technique in digital communication applications. Bandwidth efficiency, relatively good 
noise immunity, constant signal envelope, and simplicity of implementations, make this 
scheme particularly attractive for use over satellite, terrestrial radio and voiceband 
telephone channels. While system analyses pertaining to the performance of M-DPSK 
receivers abound in the literature, treatment is usually restricted to the case where 
signal transmission takes place over an additive white Gaussian noise channel. [Ref. 1] 

Since the digital information transmission is based on carriers modulated by the 
symbol waveforms, the symbol error probabilities which quantify receiver performance 
can often be computed directly as demonstrated by the many derivations of such 
probabilities found in the literature ( see [Ref 2: pp. 228-234] and [Ref 3: pp. 200-212] 
for example ). However, from the point of view of the recipient of the digital 
information, it is the Bit Error Ratio (BER) that becomes important in digital 
communication applications. While results on the BER performance of M-ary 
modulation receivers is far more limited, closed form expressions for the BER of M- 
PSK receivers have been obtained by Lee [Ref 4], with similar results independently 
derived by Tan [Ref 5]. 

This thesis is devoted to analyzing the application of direct bit detection methods 
which allow the recovery of bits in the received symbol individually, in differential 
PSK digital communication systems. For analysis purposes, 8-DPSK and 16-DPSK 
systems are considered and their performances evaluated. The particular question to 
be answered is whether for 8-DPSK and 16-DPSK modulated signals, direct bit 
detection methods provide comparable performance to conventional phase detection 
receivers, especially under conditions of improper receiver carrier phase 
synchronization. 

In this introductory chapter, background information on M-PSK communication 
systems is given. Known optimum receiver structures are presented, and Direct Bit 
Detection Receivers (DBDR) are introduced. In Chapters II and III a presentation of 
the transmitter and receiver structures respectively is given along with their operational 
properties, and the logic implemented in the DBDR for 8-DPSK and 16-DPSK 
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modulated signals is discussed in detail. Chapter IV includes the noise performance 
analyses of 8-PSK and 8-DPSK receivers, and Chapter V presents similar receiver 
analyses for 16-PSK and 16-DPSK. In both of these chapters, receiver local oscillator 
phase errors have been incorporated so as to be able to account for degrading effects 
due to improper carrier phase synchronization. The obtained mathematical expressions 
of receiver performance are evaluated on a computer and graphically presented as 
plots of bit error ratio as a function of signal to noise ratio (SNR) and phase errors. 

A. BACKGROUND 

Direct bit detection methods utilize two signal processing channels in the receiver 
in order to recover the digital information without phase angle measurements. 
Moreover, these methods exhibit optimum performance in a BER sense when carrier 
synchrOxi-iZution is completely achieved. Because of the attractiveness and simplicity of 
DBDRs, some studies have been carried out in this area, and a modified two-channel 
receiver for 8-PSK has been built and its noise performance measured by Thompson 
[Ref. 6]. Furthermore, noise performance evaluations of DBDR for 8-PSK have been 
presented by Myers and Bukofzer [Ref. 7], while similar analyses for the 16-PSK case 
can be found in Reference 8. Advantages of direct bit detection methods are well 
documented in these referenced works, in which all analyses have been done assuming 
no receiver local oscillator phase error, Gray Code bit to symbol mapping and an 
additive white Gaussian noise interference model. We present here 8-PSK, 8-DPSK, 
16-PSK and 16-DPSK communication systems under similar assumptions, except now 
the receiver has an assumed known local oscillator phase error. 

B. THEORY OF M-PSK COMMUNICATION SYSTEMS 
1. Representation of M-PSK Signals 

For M-PSK modulation with equal signal energies, a convenient 
representation [Refs. 3: p.p. 192-205] of the signal set is given by 

S;(t) = V2E7f3Cos[2rtf^t + 2;t(i-l)/M], i= 1,2,...,M , 0 <t <T^ (1.1) 

Ts 

E3= J S;2(t)dt , i=l,2,..,8 (1.2) 

0 

Tj = symbol duration 

fp = carrier frequency (Hz.) 

12 



for convenience > 1/T is assumed, 
o s 

A suitable orthonormal signal set for the representation of the above signals is 

(pj(t)= VFT COS(2jTf^t) 

(p,(t)= VW, sin(27tfj) , 0 < t < (1.3) 

so that by using trigonometric identities, Sj (t) can be expanded in terms of (Pj(t) and 
(p,(t). That is, 

Si(t) = (pi(t) V2E^ cos[ 27t(i-l)/M ] - (p 3(0 sin[ 2rt(i-l)/M ] (1.4) 

where i=l,2 ,M . 

A plot of the signal constellation for M-PSK is shown in Figure 1.1 for 
various values of M, where the coordinates of the i^^ signal vector representing the 
signal Sj(t) are 

cos[27t(i-l)/M] } ,- sinI27t(i-l)/M] } (1.5) 

For M-DPSK, the signal representations are exactly same as those already presented, 
except that the digital information is transmitted in terms of phase differences between 
consecutive signal transmissions, rather than absolute signal phases. 

2. Gray Coding for M-PSK Signals 

Before proceeding with the Gray coding for M-PSK signals, we introduce 
some definitions which will henceforth be used 

1. Symbol : Transmitted or received signal in sinusoidal form (Sj(t) where 

i= 1,2,...M) which is related to the symbol state. 

2. Symbol stale : The binary digit assignment to the symbols which always 
consists of a block of log 2 M bits in length(Sj where i= 1,2,.. .,M.) 

3. Data state : Block of log 2 M bits corresponding to the data that must be 
delivered to the intended user, and which is processed at the encoder input to 
determine the symbol to be transmitted. {D.^ w’here i= 1,2,... .,M.) 

The relationship between symbol error probability and bit error probability in 
the corresponding log 2 M bit groups depends upon the assignment of bits to the symbol 
states. The preferred bit to symbol state mapping in most cases is the Gray Code, in 
which the symbol states "which correspond to adjacent phases differ in only one bit 
position. This is because in the coherent demodulation process the most likely error 
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Figure 1.1 Signal Constellations for (a) BPSK (b) QPSK (M = 4) (c) 8-PSK (d) 16-PSK. 
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involves mistaking a correct symbol with its nearest neighbor and so such symbol 
errors translate into single bit errors. 

In Direct Bit Detection Receivers (DBDR), Gray code assignment to the 
symbol states are selected for the above given reason. In this thesis, a modified Gray 
code is used, which for comparison purposes is shown along with the conventional 



TABLE 1 

COMPARISON OF CONVENTIONAL AND MODIFIED GRAY 

CODING SCHEMES 



Svmbol 

States 


Grav 

Coding 




.Modified Gray 
Coding 




8-PSK 


16-PSK 


8-PSK 


16-PSK 


Si 


000 


0000 


111 


nil 


S2 


001 


0001 


101 


1101 


S3 


Oil 


0011 


001 


1001 


S4 


010 


0010 


on 


1011 


S5 


110 


0110 


010 


0011 


S6 


111 


0111 


000 


0001 


S- 


101 


0101 


100 


0101 


Ss 


100 


0100 


no 


0111 


S9 




1100 




0110 


Sio 




1101 




0100 


Sii 




nil 




0000 


S12 




1110 




0010 


Sl 3 




1010 




1010 


Sl 4 




1011 




1000 


Sl 5 




1001 




1100 


S16 




1000 
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Gray Code in Table I. 

Using the modified Gray Code on Table 1, the vectors representing the 
symbols with components along the (pj and (p 2 have been shown (see Figure 1.2) 
such that rj is the symbol component along the <Pj axis and r 2 is the symbol 
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component along the (P 2 axis. Therefore for M = 8 , (see Figure 1.2.a) in order to 
recover Most Significant Bit (LB) of the transmitted symbols correctly, r^ > 0 whenever 
the symbol states 5j, § 2 , Sj and Sg are transmitted and rj <0 whenever the symbol 
states S., S^, and are transmitted. The Middle Bit (MB) of the transmitted 
symbols can be recovered correctly if (r^)^ - (r 2 )^ > 0 whenever the symbol states 5j, 
S 5 and Sg are transmitted and (r^)" - (r 2 )^ < 0 whenever the symbol states § 2 , Sy 
5^ and are transmitted. The Least Significant Bit (RB) is recovered correctly if 
r 2>0 whenever the symbol states § 2 , Sy and are transmitted and r 2<0 
whenever the symbol states Sy Sy S^, and Sg are transmitted. 

Similar logic has been set up for the M= 16 case (see Figure 1.2.b) but now 
each symbol is related to 4 bits in a symbol state. Therefore for this case, 4 bits must 
be recovered from each received cy’^bol, where the bits in a symbol state are labeled 
Most Significant Bit (LB), Middle Bit Right (MBR), Middle Bit Left (MBL), and Least 
Significant Bit (RB). For example from Figure l.2.b, Sj is 0101, corresponding to 
LB = 0, MBL=1, MBR=0 and RB=1. In order to recover the LB correctly, it is 
necessary that rj > 0 for transmitted symbol states -1’ -2’ -3’ -4’ ^13’ -14’ ^15’ and 
and that r^ < 0 for transmitted symbol states Sy Sy Sj, Sg, Sg, 5jg, 5jj, and 5 j 2 * 
If r 2>0 whenever the symbol states Sj, § 2 , Sy S^, Sy Sy 5^, and 5g are transmitted 
and r 2<0 whenever the symbol states Sg, 5jg, Sj 2 , 5^^, and are 

transmitted, the RB will be recovered correctly. If (tj)^>(r 2 )^ for transmission of 
symbol states 5p S 2 , Sj, Sg, Sg, 5jg, 5jj, and and (rj)^<(r 2 )^ for transmission of 
symbol states Sy S^, Sy Sy 5 j 2 » and the MBL will be recovered 
correctly. Finally, if [ (rj)^ - (r 2 )^ - (2 rjt 2 )^ > 0 upon transmission of the symbol 

states 5p 5^ , Sy Sg, Sg, S^y ^, 3 , and 5j^, and [ (rj)^ - (r 2 )^ f - (2 x^t 2 f ^ 
transmission of the symbol states on the states S 2 , Sy Sy S^, 5jq, 5p, 5j^, and 
then the MBR will be recovered correctly. 

In subsequent chapters, the receiver logic and its performance will be 
explained in detail. Additionally, application of this methodology will be presented in 
so far as differential PSK modulation is concerned. Before doing so however, a 
summary of known results on the performance of conventional and DBDR's for M- 
PSK and M-DPSK is presented next. 

3. Known Results on the Performance of M-PSK and M-DPSK Receivers 

In this section, known results on the receiver performance for M-PSK and 
DPSK modulated signals are presented for two cases, namely conventional receivers 
and Direct Bit Detection Receivers (DBDR). 
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Figure 1.2 Symbol State Assignments for (a) 8-PSK (b) 16-PSK. 
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a. Conventional M-PSK and M-DPSK Receivers 

Coherent demodulation of M-PSK implies in principle the use of M signal 
processing channels in the receiver, as shown in Figure 1.3 (from Ha [Ref 10] ) where 
different forms of the optimum receiver structure of M-PSK signaling have been given. 
Performance analysis results for optimum receivers are weU-documented in terms of 
symbol error rate as a function of signal to noise ratio [Refs. 2,9], while more recent 
analyses have focused on the receiver probability of bit error (see [Ref 4] and 
[Ref 3: p.p. 204-207 ] for example). Symbol error probabilities can be computed from 
knowledge of the channel and signal characteristics. However for transmission of 
binary data, when analyzing system performances as a function of different levels of 
modulation, the bit error probability rather than the symbol error probability is of 
interest, as previously explained. In Prabhu [Ref 1: p. 198], the bit error prrbJiility 
has been presented only for the case of M-ary orthogonal signal sets. However, since 
the .M-PSK signal set is not orthogonal, those results are not applicable here. The 
closed form expression for the BER of M-PSK has been determined, when a Gray 
Code bit mapping is used [Ref 2: p.l98]. The results can be summarized as follows 

n 

Pb(M) = Z (k/n) P^ (1.6) 

k= I 

where 

M = Number of distinct signal waveform, 

P^ = Probability of k bits in error in a received n bit data 
block [Ref 5: pp.24-32], 
n = log 2 .M , 

Pb(M) = Bit error probability. 

A derivation of Equation 1.6 and its numerical evaluation (see Figure 1.4 ) has been 
carried out by Tan [Ref 5; p.33j. For .M-PSK signalling, receiver performance analyses 
in terms of bit error rate have also been presented by Lee [Ref 4: p.491], and plotted in 
Figure 1.5. 

b. Direct Bit Detection Receivers (DBDR's) 

DBDR structures for M-PSK have been derived for various values of .M 
[Refs. 6,7,8]. As previously indicated, this method is particularly well suited for M-PSK 
communication systems. The DBDR structures for 8-PSK (see Figure 1.6) and I6-PSK 
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(a) 




lt>) 



Figure 1.3 Block Diagram of Optimum Receiver (a) Correlator Realization 
(b) Matched Filter Realization (Ironi Ha [Ref. 10|). 
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BIR ERROR PROBABILITY LOG PIE) 




Figure 1.4 Bit Error Probabilities Versus E^Nq for 16-PSK, 8-PSK and QPSK (Ref. 5 p. 
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(see Figure l.S) have been presented and their performances have been derived in 
References 6,7,8. The BER is quantified by the probabilities of a bit error, denoted PE. 

For 8-PSK, PE has been derived by Myers and Bukofzer [Ref 7: p.l3] 

where 

PE= 1/3 { Q[ 2 VSi^ cosa ] + Q[ 2 VSNR sina ] + 

+ Q[ 2 VSNR ( cosa+ sina) ] + Q[ 2 VSNR ( cosa- sina) ] - 
- 2 Q[ 2 Vsi^ ( cosa - sina) ] . Q[ 2 VSNR ( cosa+ sina) ] } (1.7) 

where 

SNR = E^/2No=3EV2Nq 
E j = Symbol energy 

Ng= Noise pow’er spectral density level (one sided) 
a = Basic signal phase angle (see Figure 1.2). 



and 



00 

Q(x) = J [ 1/V^] exp(-u2/2) du (1.8) 

X 

The derivation of Equation 1.7 assumes an AWGN channel. Gray coded symbols and 
equal a priori bit probabilities (i.e. , 1/2). Receiver performance in BER sense is 
plotted in Figure 1.7 for various values of a as a function of SNR. 

For the 16-PSK case, probability of bit error is derived by Bukofzer 
[Ref 8] and the resulting receiver performance is plotted in Figure 1.9 for various 
values of the angles a and p, where p is shown on Figure 1.2.b. 

In both analyses, it has been assumed that the receiver's local oscillator 
does not have any phase error. In this research, similar analyses to those conducted 
in References 7 and 8 have been carried out by considering the effect of a receiver 
local oscillator phase error for 8-PSK and 16-PSK cases, as well as the noise 
performance of 8-DPSK and 16-DPSK DBDR, presented as BER as a function of 
SNR. 
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Figure 1.6 Block Diagram for 8-PSK DBDR. 
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Figure 1.8 Block Diagram for 16-PSK DBDR. 
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II. TRANSMITTER 



The generation of differential PSK signals can be performed by pre-processing 
the data bits. This chapter is devoted to anah^zing the M-DPSK transmitter structure. 
One general form of the transmitter in block diagram form for DPSK communication 
systems is shown in Figure 2.1. E.xplanation of the general elements in Figure 2.1 wUl 
be included along with the derivation of the transmitter structure for 8-DPS K and 
16-DPSK modulated signals. 




Figure 2.1 DPSK Transmitter Block Diagram. 



A. TRANSMITTER STRUCTURE FOR 8-DPSK AND 16-DPSK MODULATED 

SIGNALS 

The main task in generating the DPSK transmitter involves the design of the 
encoder, which depends on the algorithm used to relate data states to the symbol states 
prior to transmission as a sinusoidal waveform. Basically, the transmitted symbol is 
the quantized equivalent of the angle between the successively transmitted symbols. 
For reasons previously e.xplained, a Gray coding scheme is utilized in the symbol state 
to symbol mapping for the two special cases of M-DPSK signalling considered here, 
namely M = 8 and .M = 16. 
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The symbol difTerential phase relative to the data state to be transmitted is given 
on Table 2 for 8-DPSK and on Table 3 for 16-DPSK modulated signals. In both 
Tables, it is assumed that an equiphase arrangement for the transmitted symbols is 
used. 







TABLE 2 




GRAY CODE REPRESENTATION OF PHASES FOR 8-DPSK SIGNAL 

SET 


LB 


MB 


RB 


Differential 

Phase 


0 


0 


0 


0 


0 


0 


1 


a 


0 


1 


1 


la 


0 


1 


0 


3a 


1 


1 


0 


4a 


1 


1 


1 


5a 


1 


0 


1 


6a 


1 


0 


0 


la 



As mentioned earlier, the main purpose of the encoder is to generate a mapping 
algorithm whereby the phase difference of consecutive symbol outputs is related to the 
input data in accordance with the information in Table 2 and Table 3. In order to 
achieve successive signal transmissions the following algorithm is implemented in the 
transmitter 

1. Store the present data state (Dj, i= 1,2...,M) and the previous symbol state 
(Si,i=I,2,3..,M). 

2. From the Table 2 or Table 3, find the corresponding differential phase value 
associated with the present data state as Na, where N is an integer v\’ith values 
ofO<N<M-I. 

3. Generate a symbol state, such that the corresponding symbol will have a phase 
angle Na degrees larger ( in additive sense, modulo 2;t) than the previously 
transmitted symbol. 
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TABLE 3 

GILAY CODE REPRESENTATION OF PHASES FOR 16-DPSK SIGNAL 

SET 



LB MBL MBR RB 

0 0 0 0 

0 0 0 1 

0 0 11 
0 0 10 
0 110 
0 111 
0 10 1 

0 10 0 

110 0 
110 1 
1111 
1110 
10 10 
10 11 
10 0 1 

10 0 0 



DifTerential 

Phase 

0 

a 

2a 

3a 

4a 

5a 

6a 

la 

8a 

9a 

10a 

11a 

12a 

13a 

14a 

15a 



For example, assuming the previous symbol state is (0101) and present data state is 
Dg (0110) which corresponds to the differential phase value of 4a (see Table 3.)- From 
Figure 1.2.b, proceed 4a degrres counterclockwise from S^, resulting in the symbol 
state is (1111) will be the sy’mbol output from the encoder. 

This algorithm can be shown in matrix structure, in which the row entries and 
the column entries of the matrix are the symbol states and data states respectively, and 
the elements of the matrix are the symbol state which correspond to the appropriate 
phase difference in the transmitted symbol (see Figure 2.3 and Figure 2.2). For 
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example, under the same assumptions given in the previous paragraph, by using 
(0110) as a row entry and 5- as a column entry in the matrix in Figure 2.3, present 
symbol output from encoder is found as (1111) at the intersection of the row and 
column of the matrix. 



DATA 


^6 

000 


^7 

100 


STMBOL STATE 
^8 ^1 ^2 ^3 

110 111 101 001 


^4 

on 


^5 

010 


STATE 
D^- 000 


000 


100 


110 


111 


101 


001 


on 


010 




100 


100 


110 


111 


101 


001 


oil 


010 


000 




110 


110 


111 


101 


001 


0x1 


vJlO 


000 


100 


111 


111 


101 


001 


oil 


010 


000 


100 


no 


dJ- 


101 


101 


001 


oil 


010 


000 


100 


no 


111 




001 


001 


oil 


010 


000 


100 


no 


111 


101 




oil 


oil 


010 


000 


100 


110 


111 


101 


001 


?5- 


010 


010 


000 


100 


110 


111 


101 


001 


on 



Figure 2.2 Encoder Matrix for 8-DPSK Signals. 

The hardware implementation of the encoder structure can be accomplished in 
several forms. A practical approach to building the encoder can be achieved by the use 
of EPRO.Ms (Erasable Programable Read Only Memory)- Once the truthe table the 
programming of the EPROM is a trivial task to accomplish. Demultiplexers can also 
be used to build the encoder circuitry; however, because such a device would be hard 
wired, it would not have the reprogrammability of the EPRO.M implementation. 

In this research, the 8-DPSK and 16-PSK communication systems are evaluated 
using a computer. Appendix A includes the source file which is used to generate the 
possible outputs of an .VI-PSK and an .VI-DPSK communication system. 

The carrier modulator uses the encoder outputs (symbol states) as an input, and 
generates the symbols (signals) suitable for transmission. The signal transmission 
involves sinusoidal voltages whose phase is related to the symbol states (encoder 
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Figure 2.3 Encoder Matrix for 16-DPSK Signals. 
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outputs) and of carrier frequency suitable for the transmission channel. The signal 
generation can be achieved via a multiplexer, whose structure is shown in Figure 2.4. 
For this kind of implementation. (M-1) phase shift elements are needed with set phase 
shift values. Another approach would involve an implementation using a D/A 
converter along with a phase modulator which can be obtained with a VCO (see 
[Ref. 11: pp. 269-27 IJ. 




Figure 2.4 Carrier .Modulator Structure for M-DPSK Signaling. 



B. COMPUTER GENERATION OF THE TRANSMITFER OUTPUTS 

For computer evaluation of the M-ar\' DPSK communication system , a 
FORTRAN program called BELIZ has been written (see Appendix A). The program 
generates all possible data state and symbol state combinations for M-DPSK, and 
produces the corresponding symbol states outputs. The program outputs are given in 
Appendix B in tabular form along with the transmitter matrix of Figure 2.3 and Figure 
2.2, specified in decimal form. 
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III. RECEIVER 



A. RECEIVER STRUCTURE FOR 8-DPSK AND 16-DPSK MODULATED 

SIGNALS 

The imroductor>’ chapter presents an over\’iew of the structure and performance 
of optimum receivers for M-ary signalling. The focus of this chapter is on DBDR's. 
.As mentioned in Chapter II, the symbol state generation based on the data state and 
the previous symbol state is the most important task performed by the encoder. 
Therefore in the receiver's decoder, a similar algorithm must be performed in reversed 
sense, that is, from two successive symbol states, the corresponding data state must be 
estimated. .A generic form of the receiver in block diagram structure is shown in Figure 
3.1. 




Figure 3.1 General receiver block diagram for DPSK signals. 

The algorithm for recovering the transmitted data is as follows 

1. Recover the symbol state from the received symbol using direct bit detection 
techniques. 

2. Store the past symbol state 5j,i= 1,2.. .,M. 

3. Compare the present symbol state with the past symbol state to determine the 
incremental " phase angle " between the two. 

4. From Table 2 or 3 find the corresponding data state that gave rise 

to the "phase angle" between the two consecutive symbol states. 
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For example, assuming the past symbol state is (S^) 0001 and the present symbol state 
is Sj (1001) from Figure 1.2. b, the angle between and is 3a. From Table 3., the 
corresponding Gray Code representation for 3a is found to be (0010) , therefore the 
recovered data state will be Dp. 

The decoder implementation can be depicted as a matrix structure in a manner 
similar to that shown for the transmitter encoder in Chapter II. For 8-DPSK and 
I6-DPSK modulated signals, the decoder matrices are given in Figure 3.2 and Figure 
3.3 respectively. For example, assuming the same conditions as in the above 
paragraph, by using and as a column and row entries for the matrix given in 
Figure 3.3, the recovered data state from the decoder is determined as Dp (0010). 









SYMBOL 


STATE 








-^6 


^7 


^8 


^1 


^2 


§3 


?4 


is 


SYMBOL 


000 


001 


Oil 


010 


110 


111 


101 


100 


STATE 


















S, - 000 


000 


001 


Oil 


010 


110 


111 


101 


100 


S, - 100 


100 


000 


001 


oil 


010 


110 


111 


101 


si - 101 


101 


100 


000 


001 


oil 


010 


110 


111 


i? - 111 


111 


101 


100 


000 


001 


oil 


010 


110 


s: - 101 


101 


111 


101 


100 


000 


001 


oil 


010 


- 010 


010 


101 


111 


101 


100 


000 


001 


oil 


K - 0^^ 


oil 


010 


101 


111 


101 


100 


000 


001 


- 001 


001 


oil 


010 


101 


111 


101 


100 


000 



Figure 3.2 Decoder matrix for 8-DPSK signals. 



B. COMPUTER GEiNERATION OF RECEIVER OUTPUTS 

Appendix C includes the results of the computer evaluation of the .M-ar\^ DPSK 
receiver. In addition to generating the decoder matrix in decimal form, the output for 
this program includes all possible combinations of the decoder outputs. 
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Figure 3.3 Decoder matrbc for I6-DPSK signals. 
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C. RECEIVER LOCAL OSCILLATOR PHASE ERRORS 



The introductory chapter summarizes the studies done on M-PSK DBDR 
structures, assuming that the receiver local oscillator has no phase error. In this 
research, a fixed receiver phase error has been assumed in the performance analyses 
conducted, and its effect on receiver BER is presented in Chapter IV and Chapter V. 
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IV. PERFORMANCE ANALYSES OF 8-PSK AND 8-DPSK DBDR'S WITH 
LOCAL OSCILLATOR PHASE ERROR 

In previous chapters, DBDR's are introduced and their operational properties are 
explained. In this chapter, performance analyses for 8-PSK and 8-DPSK will be 
carried out in detail. 

A. 8-PSK DBDR PERFORMANCE ANALYSIS IN THE PRESENCE OF 
RECEIVER LOCAL OSCILLATOR PHASE ERROR 

The possible transmitted signals which are represented as vectors in Figure 1.2.a , 
can be mathematically expressed as 

Sj(t) = V2E7f; sin [ 27tfot + 0j(t) ] , i= 1,2,3...,8 , 0< t ^T^ (4.1) 

where 



E^= J Sj2(t)dt , i=l,2,..,8 (4.2) 

0 

T = svmbol duration 
s 

fjj = carrier frequency (Hz.) 

6j(t) indicates the phase of the transmitted PSK signal, which for the 8-PSK case 
becomes 



e-.(t) = 



a + (i-1) 71/4 i= 1, 3,5,7 
i7i/4 - a i = 2,4,6,8 



0<t<Tj 



(4.3) 



The angle a is variable as shown in Figure 1.2. a; however, the constraint 
0 a ^ 7t/4 must be observed. Equation 4. 1 can be written by using a suitable 
onhonormal set of functions as 



Sj(t)= <Pi(0 VEs cos0j(t) + (p 2 (t) VEj sin0;(t) 



(4.4) 
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i=l,2,...,S : 0<t<T^ 

where (pj(t) and cp 2 (t) are the orthonormal function set which by applying the 
Gram-Schmitt orthonormalization procedure 
[Ref 3: p.l98] 
can be shown to be 



(pj(t) = V2/T sin 27tfgt 

(P2(t) = VItTcos 2;tfQt , 0< t <T^ (4.5) 

The received signal depends upon the channel model and the transmitted signal 
as specified by Equation 4.1. In this research, an Additive White Gaussian Noise 
(AWGN) interference model has been assumed so that the received signal can be 
written as 

r(t) = Sj(t) + n(t) (4.6) 

where n(t) is a sample function of a white Gaussian noise process having a constant 
(one sided) power spectral density function equal to Nq watts/Hz. The receiver 
processes r(t) over an interval of duration T^ seconds and a decision is made at the end 
of the interval. This procedure continues in a repetitive fashion as time progresses. 

Assuming all signals are equally likely to be transmitted, the probability of a bit 
being correct for M-PSK signal transmission is written as 

M 

Pr{ bit correct } = l/M ^ Pr( bit correct / Sj(t) transmitted } (4.7) 

i= 1 

The calculations for the conditional probabilities (Pr{ bit correct / Sj(t) transmitted }, 
i= 1,2,.... M ) have been determined by Myers and Bukofzer [Ref 7: pp. 6-12], for the 
8-PSK case with no receiver local oscillator phase error. In the sequel, the same 
analysis methods are utilized to obtain the receiver BER in the presence of a local 
oscillator phase error. A local oscillator phase error is modeled by the angle £ in the 
unmodulated receiver carriers cpj(t) and ^ 2 ^^) Figure 1.6. These are denoted (pj'(0 
and <p 2 '( 0 » namely 
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9i’(t) = V2/T sin ( 2;:fgt + e ) 

(p,'(t) = V27T cos ( 2;rfQt + c ), 0 < t < 



(4.8) 



.Assuming signal Sj(t) was transmitted, this produces outputs j from the receiver's 
correlators, where 



^k,i = J [ S;(t) + n(t) ] (p^'(t) dt , i= 1,2,...8 , k= 1,2 
0 



(4.9) 



For example when Sj(t) is transmitted, from Equation 4.4, 6.(t) becomes just a and 



rj 2 = VEj (®’^) 
rj 2 = VEj sin (a-e) + U 2 



(4.10) 



where 



T 

S 

n^ = J n(t) (pj^'(t) dt , k= 1,2 
0 



(4.11) 



Observing that 



T, T 

S S 

E [nj n 2 l = E [ J n(t) (pj'(t) dt x J n(t) dt ]= 0 
0 0 



(4.12) 



This means that nj and n 2 are uncorrelated random variables. Since nj and n 2 are 
Gaussian random variables (having identical variances) [Ref 7: p.7], they are also 
statistically independent. 

The goal of the receiver of Figure 1.6 is to correctly recover the LB, the MB, and 
the RB, transmitted via the signal Sj(t). By introducing the more compact notation 



Pr{ bit correct / Sj(t) } = Pr{ bit correct / Sj(t) transmitted } , i = 1,2. .,8 (4.13) 
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the probabilities of such correct bit recovery are given in general by 



Pr{ LB correct / Sj(t) } = i 



Pr{ MB correct / Sj(t) } = < 

\ 

/ 

Pr{ RB correct / Sj(t) } = < 



Pr{ri^j>0/Sj(t)} i= 1,2, 8, 7 

Pr{ rj j<0 / Sj(t) } i= 3,4,5, 6 
Pr{ri.2-r2,.2>0/Si(t)} i= 1, 4,5,8 

2,3,6,7 

Pr{r2^i>0/S;(t)} i= 1,2, 3,4 
Pr{ r 2 ,j < 0 / Sj(t) } i= 5,6,7,8 



(4.14) 



(4.15) 



M 16 ) 



From Equations 4.8-4.16, we can write 



Pr{ LB correct / Sj(t) } = Pr( r^ j > 0 / Sj(t) } 

= Pr{ nj > - cos(a-c) / Sj(t) } 

= 1 - Q [ V2 E^/Nq cos(a-c) ] (4.17) 

where 

00 

Q(x) = J [ 1/V2jt] exp(-u^/2) du (4.18) 

X 



Similarly 

Pr{ RB correct / Sj(t) } = Pr( > 0 / Sj(t) } 

= Pr{ n^> - sin(a-e) / Sj(t) } 

= 1 - Q [ V2 Ej/Nq sin(a-e) ] (4.19) 

In order to evaluate the conditional probability of the MB being correct, the following 
steps are required. From Equation 4.15 



Pr{ MB correct / Sj(t) } 



M,i ‘2,i - 

Pr{(ri,.-r 2 ,.).(ri.+ r^.) >0 / s. 



(4.20) 
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By defining variables j and y^ . as 



>’l,i = + •‘2,1 

>':.i " ^i,i * h,i 



(4.21) 



then from Equation 4.9, when Sj(t) is transmitted, 

^ 1,1 ~ t ''' ^^2 

>’2 2 = V Ej [ cos (a-£) - sin (a-e) ]+ Hj - n 2 



(4.22) 

(4.23) 



Clearly y^ 2 and ^ are Gaussian random variables with statistical means given by 



E{ yj [ cos (a-£) + sin (a-£) ] 

E{ >’ 21 ) = n/^ [ cos (a-£) - sin (a-£) ] 



(4.24) 

(4.25) 



and variances given by 



Var{yi,i} = Var{y2^^} =N, 



(4.26) 



By observing that 



E { (yi_i- E [ y^^J) (y2^j- E [ y2^ J) } = E { (n^ + n2) (n^- n2) } 

= E { (n^)2-(n2)2 } = 0 (4.27) 

where n^ and n 2 are defined in Equation 4.11. Since yj ^ and J 2 j are uncorrelated 
Gaussian random variables, they are also statistically independent. Therefore, 
Equation 4.20 becomes 

Pr{ MB correct / Sj(t) } = Pr{ y^ ^ y^^^ >0 / s^(t) } 

= Pr{ [ yi,i > 0 / s^(t) ] and [ y 2 ^^ > 0 / s^(t) ] ) + 

+ Pr{ [ y^ j < 0 / Sj(t) ] and [ y 2 ^j < 0 / s^(t) ] } + 

= Pr{ yj^j > 0 / Sj(t) ) Pr{ y^ ^ > 0 / s^(t) } + 

+ Pr{ y^’i < 0 / Sj(t) } Pr{ y 2 ’i < 0 / s^(t) } (4.28) 
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From Equations 4.22 - 4.26, we have 

(4.29) 

(4.30) 



Pr(vi^i>0/si(t)} = Q[E{yi^j}/ VN q] 

Pr( 724 >0 / } = Q [ £{724} / ] 



so that Equation 4.28 becomes 

Pr( MB correct / s^(t) } = 1 - Q{ V E^/Nq [ cos(a-£) + sin(a-c] } - 

- Q{ t cos(a-£) - sin(a-£] } 

- 2 Q( V £ 5/^0 f cos(a-£) - sin(a-£] } . 

. Q{ V t cos(a-£) 4- sin(a-£] } (4.31) 

Assuming now transmission of s^(t), a similar analysis procedure is performed 
but now Equations 4.14-4.16 become 



02 (t) = 7 t /2 - a 


(4.32) 


rj 2 = sin (a+ £) 4- nj 

’*2 2“ cos(a4-£) 4- n 2 


(4.33) 



where n^ and n 2 are defined by Equation 4.11. From Equations 4.14 and 4.33 



Pr{ LB correct / S 2 (t) } = Pr{ 2 ^^ / ^ 2^0 ) 

= Pr{ n 2 > - sin(a + £) / S 2 (t) } 

= 1 - Q [ VTE^Nq sin(a + £) ] (4.34) 

and 



Pr( RB correct / S 2 (t) } = Pr{ ^ 22 ^^ ^ ^ 

= Pr{ n 2 > - cos(a + £) / S 2 (t) } 

= 1 - Q [ J2 E 3 /N 0 cos(a + £) ] (4.35) 

For calculations involving the MB, one additional random variable needs to be defined, 
namely 



^ 3 , 



''2,1 ■ ’’l,i 



i= 1,2,....8 



(4.36) 
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Assuming transmission of s,(t), Equation 4.36 becomes 
'3.2 ^2,2 ” ^1,2 

= V Ej COS (a + e) + n2 - V E^ sin (a + £)-n^ 



(4.37) 



where 

E (vj ,} = cos (a+£) - sin (a + £) = E (4.38) 

Since 

Pr( MB correct / s.,(t) } = Pr{ (rj ^)^ - ^ ^ (^-39) 

— Pr{ (fjj ■ ^ 2 , 2 ) ^^ 1,1 ^ 2 , 2 ^*^^ ^ ^ 

by using Equations 4.21 and 4.36 

Pr( MB correct / s,(t) } = Pr{ y^^^ 3^3,2*^ ^ ^ 

Furthermore 

E { (>’1^2* E ^ 3 'i, 2}) (>’3,2* E { 73,2)) ) = E { (1^1 + n^) (n^- n^) } 

= E ( (ni)2-(n2)^ } = 0 (4.41) 

so that y^ 2 y^ ^ statistically independent Gaussian random variables, and 
therefore Equation 4.40 becomes 

Pr( MB correct / S 2 (t) } = Pr{ y 3^2 ^ 1,2 ^ ^ ^ 

= Pr{ I yi, 2 > 0 ! Sj(t) ] and [ 0 / S 2 (t) ] } + 

+ Pr{ ( yi^2 ^ ^ / ^2^^) 1 I >^3,2 ^ 1 ^ ^ 

= Pr{ yi^2^ ^ ^ ^2^^) ^ >^3,2^ ^ I ® 2 (t) } + 

+ Pr{ yi,2*^^ / ^2(0 } Pr{ y^,2^^ t ^2(0 } ( 4 - 42 ) 

By using a method similar to that used in the derivation of Equations 4.29 and 4.30, 
for yj 2 and y^ ^ » Equation 4.42 becomes 
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Pr{ MB correct / s,(t) } = 1 - Q{ V Ej/Nq [ cos(a + £) + sin(a + £] } - 

- Q{ V Ej/Nq [ cos(a + £) - sin(a + £] } 

- 2 Q{ V t cos(a + £) + sin(a+£] } . 

• Q{ V Ej/Nq [ cos(a + £) + sin(a + £] } (4.43) 

Analyses can be carried out for probability of correct reception assuming every 
possible signal transmission. This is however not necessary because of inherent 
symmetries in the problem . In conditions where symmetries are not present, previous 
methodologies can be used to obtain results as follows. Assuming transmission of 
83 ( 0 , Equations 4.3 and 4.9 become 



03 (t) = Jt /2 + a 


(4.44) 


tj 3 = cos (a-£+7t/2) + nj 


(4.45) 


= - VEj sin («-£) + n^ 




r^ 3 = VEj sin (a-£ + ;t/ 2 ) + 


(4.46) 


= - VEj cos (a-£) + ri2 





so that from Equation 4.14 

Pr{ LB correct / Sj(t) } = Pr{ Tj 3 < 0 / Sj(t) } 

= Pr{ nj > - sin (a-£) / S 3 (t) } 

= 1 - Q[ J2 E^/Nq sin(a-£) ] (4.47) 

and from Equation 4.16 

Pr{ RB correct / S 3 (t) } = Pr{ r^ 3 < 0 / S 3 (t) } 

= Pr{ n 2 > - VE^ cos (a-£) / 83 c t) } 

= 1 - Q[ VTeTNo cos(a + £) ] (4.48) 

Assuming transmission of s^(t), Equations 4.3 and 4.4 become 

= K- a (4.49) 

= V^s (Jt-a-£) + (4.50) 

= - cos (a + £) + nj 
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(4.51) 



4 ~ (?t-a-£) + n, 

= - sin (a + €) + n2 

so that from Equation 4.14 

Pr{ LB correct / s^Ct) } = Pr{ < 0 / s^(t ) } 

= Pr{ n^ > - sin (a + c) / s_^(t) } 

= 1 - Q[ V2“e7No cos (a + c) ] (4.52) 



and from Equation 4.16 

Pr{ RB correct / s_j(t) } = Pr{ r 2 4 > 0 / S 4 (t) } 

= Pr{ n.,> - (®'^^) ' ^ 4 (t) } 

= 1 - Q( VTeTNq sin (a + £) ] (4.53) 



Also from Equation 4.15 

Pr{ MB correct / S4(t) } = Pr{ (rj^4)^ - (r2^4)^> 0 / S4(t) } (4.54) 

= Pr{ (rj ^4 - r2^4) (r ^^4 + ^ ^ 

Using Equations 4.21 and 4.36 yields 
>'1,4 = ^1,4 + ^4 

'2,4^ ^1,4 '^2,4 

and using a method similar to that used in the derivation of Equation 4.27 results in 

Pr{ MB correct / S4(t) } = Pr{ yj^4 y2^4 >0 / $4(0 } 

= Pr{ [yi,4>0/ 84(0 ] and [y2^4>0/ S4(t) ] } + 

+ Pr( [ yi,4<0 / 54(1) ] and [ y2,4<0/ S4(t) ] } + 

= Pr{ yj^4> 0 / S4(t) } Pr{ y2^4> 0 / S4(t) } + 

+ Pr{ yj’4 < 0 / S4(t) } Pr{ y2*4 < 0 / S4(t) } (4.57) 



(4.55) 

(4.56) 



which becomes 
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Pr( MB correct / s^(t) } = 1 - Q { V2 E^/Ng [ cos(a + e) - sin(a + e) ] } - 

* Q { V2 [ cos(ct + £) + sin(cH- £)•,]} + 

+ 2 Q { V2 E^/Ng [ cos(a + £) - sin(a + £) ] } . 

. Q { Jl Ej/Nq [ cos(a + £) + sin(a + £) ]} (4.58) 

For transmission of S5(t), the pertinent equations are 



65(0 = 71+ a 


(4.59) 


tj 5 = V^s cos (7T + a-£) + Uj 
= - cos (a-£) + Uj 


(4.60) 


’■25” (7t + a-£) + U 2 

= - -v^s '■’n (a-£) + U 2 


(4.61) 



so that 



Pr{ MB correct / S5(t) } 



= Pr{(r^ 5)2-(r2^5)2> O/SjCt)} 

= Pr{ (tj 5 - r2^5) (tj ^5 4- r 2 5>>0 / SjCt) } 



(4.62) 



which becomes 

Pr{ MB correct / Sj(t) } = Pr{ yj ^ 72 5 > 0 / S5(t) } 

= Pr{ [ yi,5> 0 / Sj(t) ] and [ 72^5> 0 / S5(t) ] } + 

+ Pr{ [ Yl^S < 0 / 85(0 ] and [ 72 5 < 0 / S5(t) ] } + 

= Pr{ yi, 5 > 0 / S 5 (t) } Pr{ Y 2 ^ 5 > 0 / S 5 (t) } -I- 

+ Pr{ 7 i’ 5<0 / Sj(t) } Pr( y2’s<0 / S 5 (t) } (4.57) 

by observing that y^ ^ and 72 5 are statistically independent random variables making 
it possible to apply the similar procedure shown above. The result becomes 



Pr( MB correct / Sj(t) } = 1 - Q { E^/Nq [ cos(a-£) + sin(a-£) ] } - 

- Q { V2 Ej/Nq [ cos(a-£) - sin(a-£) ]} + 

+ 2 Q { V2 Eg/Ng [ cos(a-£) + sin(a-£) ] } . 

. Q { V2 E^/Nq [ cos(a-£) - sin(a-£) ]} (4.64) 
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As indicated before, because of the inherent symmetry in the problem, the 
following equations hold for conditional LB correct probabilities (see Figure 1.2.b ) 



Pr{ LB correct / Sj(t) } = Pr[ LB correct / Sj(t) } 


(4.65) 


Pr{ LB correct / s,(t) } = Pr{ LB correct / s^(t) } 


(4.66) 


Pr{ LB correct / Sj(t) } = Pr{ LB correct / } 


(4.67) 


Pr{ LB correct / s^(t) ) = Pr{ LB correct / Sg(t) } 


(4.68) 


and similarly for the RB 




Pr{ RB correct / s^(t) } = Pr{ RB correct / S 5 (t) } 


(4.69) 


Pr{ RB correct / S 2 (t) } = Pr{ RB correct / s^(t) } 


(4.70) 


Pr{ RB correct / S 2 (t) } = Pr{ RB correct / S 2 (t) } 


(4.71) 


Pr{ RB correct / s^(t) } = Pr{ RB correct / Sg(t) } 


(4.72) 


and finally for the MB 




Pr{ MB correct / s^(t) } = Pr{ MB correct / Sj(t) } 


(4.75) 


Pr{ MB correct / S 2 (t) } = Pr{ MB correct / Sg(t) } 


(4.76) 


Pr{ MB correct / s^(t) } = Pr{ MB correct / s^(t) } 


(4.77) 


Pr( MB correct / S 5 (t) } = Pr{ MB correct / Sy(t) } 


(4.78) 



By using now Equation 4.7 for the unconditional probability of the LB, MB and 
RB being correct, the above results yield 



Pr{ RB correct } = Pr{ LB correct } 

= 1 - 1/4 [ Q[ V2E;7No cos(a-£) ] + 

+ Q[ V2 Ej/Nq sin(a + £) ] + 

+ Q[ V2 Ej/Nq cos(a + £) ] 

+ Q[ vTeJNo sin(a-£) ] ] (4.79) 

and 

Pr{ MB correct } = 1 - 1/2 [ Q { E^/Nq [ cos(a-£) + sin(a-£) ] } + 

+ Q { >/2 Ej/Nq [ cos(oi + £) + sin(<z + c) ]} + 
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(4.80) 



+ Q { V2 Ej/Nq [ cos(a-E) - sin(a-£) ] } + 

Q { V2 Ej/Nq [ cos(a + £) - sin(a + £) ]} ] + 

+ Q { V2 Ej/Ng [ cos(a-£) + sin(a-£) ] } . 

. Q { V2 Ej/Ng [ cos(a + £) + sin(a + £) ]} + 

+ Q { V2 Ej/Nq [ cos(a-£) - sin(a-£) ] } . 

. Q { V2 e}Nq [ cos(a + £) - sin(a + £) ]} ] 

For 8-PSK,the Signal to Noise Ratio (SNR) has been defined by Equation 1.6, so 
that using Equations 4.79 - 4.80, the BER denoted PE can be expressed as a function 
of SNR, a and £ as follows 

PE= 1/3 [ Pr{ LB correct } + Pr{ MB correct } +Pr{ RB correct } ] 

= 1 - 1/6 { Q[ 2 VSNR cos(a-£) ] + Q[ 2 Vs!^ sin(a + £) ] + 

+ Q[ 2 V SNR cos(ot + £) ] + Q[ 2 -J SNR sin(a-£) ] + 

+ Q { 2 V SNR [ cos(a-£) + sin(a-£) ] } + 

+ Q { 2 V SNR [ cos(tf + £) + sin(ct + £) ] } + 

+ Q { 2 J SNR [ cos(a-£) - sin(a-£) ] } + 

+ Q { 2 V SNR [ cos(ct + £) - sin(ct + £) ]}}■!■ 

+ Q { 2 V SNR [ cos(a-£) + sin(a-£) ] }. 

. Q { 2 J SNR [ cos(a + £) + sin(a + £) ] } + 

+ Q { 2 J SNR [ cos(a-£) - sin(a-£) ] }. 

. Q { 2 V~S!^ [ cos(a + £) - sin(a + £) ] } (4.81) 

From Equation 4.81, the PE can be computed for various values of SNR, a, and 
£, with the results used to optimize the DBDR performance. For £= 0° (i. e. , no 

receiver local oscillator phase error), values of PE vs SNR have been plotted (see 
Figure 4.1) showing complete agreement with the results obtained by Myers and 
Bukofzer [Ref 7; p. 11] Figure 4.1 includes also the PE for various values of £, as a 
function of SNR. 

B. 8-DPSK DBDR PERFORMANCE ANALYSIS IN THE PRESENCE OF 
RECEIVER LOCAL OSCILLATOR PHASE ERROR 

The advantages and disadvantages of differential encoding and decoding have 
been previously explained in the introductory chapter. In this section, performance of 
the 8-DPSK DBDR is presented. As explained earlier the 8-DPSK decoder utilizes the 
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LB, MB and RB. which are recovered via direct bit detection methods, and performs an 
inverse mapping from which the original data stream is delivered to the intended user 
Before proceeding with the performance calculations, the following notation will 
be introduced. The probabilities of a bit being correct or being in error w'ill be denoted 

Ck = Pr{ k correct } , Ek = Pr{ k error ) (4.82) 



respectively, where k= 1 for LB, k= 2 for MB, k= 3 for RB. By introducing the 
notation given by Equation 4.82, the analysis for the performance of the 8-DPSK 
DBDR can be expressed in more compact form. For example E2 means the 
probability that the MB is in error. 

In order to obtain the probability of a bit being correct (with differential 
decoding) we must first evaluate 

Pr{ RB correct } = 1/2 Pr{ RB correct / RB= 1 } + 

+ 1/2 Pr( RB correct / RB= 0 ) (4.83) 

where by further conditioning, we equivalently have 

1 1 

Pr{ RB correct/RB = m ) = Y, X correct / RB = m,MB = i,LB = j }. 

i=0 j = 0 

. Pr{ MB = i, LB = j } m=0,l (4.84) 

Due to the assumption of equally likely probabilities of having O's and I's in the data, 
Pr{ MB = i, LB=i } = 1/4 i,j= 0,1 (4.8o) 

so that Equation 4.83 becomes 

1 1 

Pr{ RB correct } = 1/8 Y Pr{ correct / RB = 0,MB = i,LB = j } + 

i = 0 j = 0 
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(4.86) 



1 1 

+ 1/8 y X Pfi correct / RB= l,MB = i,LB = j } 
i = 0 j = 0 

The DBDRs recover from the received signal, each bit in the data state (for PSK 
modulation) or each bit in the symbol state (for DPSK modulation) independently. 
This property of the DBDR is an advantage for BER analyses. Because certain 
consecutive bit errors in the past and presently recovered symbol states do not 
necessarily result in errors in the recovered data states, these bit errors in successive 
symbol states must be considered. That is, even though symbol errors may occur, the 
delivered data may stillbe correct. Therefore, these type errors have to be accounted 
for as resulting in a correct reception and considered in the computation of the 
probability of a bit being correct. Figure 4.2 shows one example of the situation being 
discussed. Here the data state (0,0,0) is assumed to occur and the symbol state (0,0,0) 
is assumed to have been previously transmitted. This means that according to the 
encoder matrix shown in Chapter II, the symbol state (0,0,0) must now be generated 
for transmission. Figure 4.2 further shows the receptions of the symbol states which 
even though, most contain errors, produce recovered data states in which the RB is 
correctly delivered. Observe that there are 32 such combinations which yield a correct 
0 for the RB. Therefore notation has to be introduced in order to account for error 
patterns of the above type, namely Pr{ RB correct, (k,l,m) / RB = 0,MB = i,LB = j }, 
where k,l,m= 0,1. The meaning of this expression is the probability of the RB 
correctly recovered given that the data state was (j,i,0) and the delivered data state is 
(k,l,m) (where m=0 in this case, otherwise the RB could not be correct due to the 
transmission of the RB = 0). Therefore 

Pr{ RB correct / RB = 0,MB = 0,LB = 0 } = 

1 1 

= Z Z RB correct,(k,l,0) / RB = 0.MB = 0,LB = 0 ) (4.87) 

k=0 1=0 



Since each element of the symbol state is independently recovered, we see from Figure 
4.2 that, for example 
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miHmmnuimmimufimimimimimimimHmimmHMmumim 

CORRECT RECEPTION COMBINATIONS IN PRESENCE OF ERROR 



INPUT PREVIOUS PRESENT PREVIOUS PRESENT 

TRANSMISSION TRANSMISSION RECEPTION RECEPTION 



000 000 000 000 000 

0 0 1 0 0 1 

0 10 0 10 

oil oil 

10 0 10 0 

10 1 10 1 

110 110 

111 111 

0 0 0 111 

0 0 1 0 0 0 

0 10 110 

oil 100 

10 0 10 1 

10 1 0 10 

1 1 0 0 0 1 

111 oil 

0 0 0 0 1 0 

0 0 1 10 1 

0 10 oil 

oil 001 

1 0 0 0 0 0 

10 1 111 

110 10 0 

111 110 

0 0 0 1 0 1 

0 0 1 10 0 

0 10 111 

oil 110 

1 0 0 0 0 1 

10 1 0 0 0 

110 oil 

111 010 



OUTPUT 



0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 10 
0 10 
0 10 
0 10 
0 10 
0 10 
0 10 
0 10 
10 0 
10 0 
10 0 
10 0 
10 0 
10 0 
10 0 
10 0 
110 
110 
110 
110 
110 
110 
110 
110 



Figure 4.2 Example of Error Patterns for 8-DPSK System Where Data State, 
Previous Transmission, and Present Transmission are 000 with RB Correctly Recovered. 
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Pr{ RB correct, (0,0,0) / RB = 0,MB = 0,LB = 0 } = Cl^ C2^ C3^ + 

+ Cl" C2^ E3-+ Cl" E2^ C3^+ Cl^ E2^ E3^+ El" C2^ C3" + 

+ El" C2^ E3^+ El^ E2^ C3^+ El^ E2^ E3^ (4.88) 

From Equations 4.79 and 4.82 - 4.86, it is clear that Cl= C3 and El= E3. 
Therefore, simplification of Equation 4.88 is possible, thus yielding 

Pr{ RB correct, (0,0,0) / RB = 0,MB = 0,LB = 0 } = Cl^ C2^+ 2 Cl^ C2^ El^ + 

+ Cl-^ E2^+ 2 Cl^ El^ E2^+ El'^ C2^+ El^ E2^ (4.89) 

A similar expression must be obtained for other combinations of values of k and 1 as 
well as the condition that the RB is 1. This produces a total of 64 possibilities of 
received symbol states that result in a correct RB. A computer generated table was 
used ,in order to identify the equivalent error patterns for other cases involving MB= 1 
and LB= 1. The appropriate listing are shown in Appendix C. Similar analyses must 
be carried out in order to determine the unconditional probabilities Pr{MB correct} as 
well as Pr{LB correct). Such analyses demonstrate that due to the inherent symmetries 

Pr( LB correct } = Pr{ MB correct } = Pr{ RB correct } = PE (4.90) 

so that using Equation 4.86 along with similar forms obtained for k = 1 and 1= 1 yield 

PE = 1 - 1/2 [ 2C1^ C2^ + 9 Cl^ C2^ El^+ Cl*^ E2^ + 3 Cl"^ E2 C2 + 

4 - 10 Cl^ El^ E2^4- 2 El'* C2^ + El"* E2^ + 17 Cl^ C2 El^ E2 4 - 

4 - Cl^ C2^ El 4- 3 Cl^ El E2^ 4 - 2 Cl Ep C2^ 4 - 2 Cl El^ E2^ 4 - 

4 - 3 C2 El"* E2 4- 4 Cl C2 El^ E2 4 - 4 Cl^ C2 El E2 ] (4.91) 

In Figures 4.4 and 4.5, PE (given by Equation 4.91 ) vs SNR have been plotted for 
various values of a and £ respectively. 

C. NUMERICAL RESULTS AND DISCUSSION OF RESULTS FOR 8-PSK 
AND 8-DPSK DBDR'S WITH RECEIVER LOCAL OSCILLATOR PHASE 
ERROR. 

Previous sections of this chapter have focused on the performance of 8-PSK and 
8-DPSK DBDRs operating in the presence of a receiver local oscillator phase error. 
Without any phase error, we have shown the performance plots of 8-PSK and 8-DPSK 
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Figure 4.3 PE vs SNR Comparison Plot for S-PSK and S-DPSK DBDR 

(a= 22.5® , £= 0®). 
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Figure 4.5 PE Versus SNR of 8-DPSK DBDR for Various Values of e (a= 22.5”). 
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DBDRs on Figure 4.3, thus demonstrating the performance penalty of the latter 
receiver. While usage of 8-DPSK provides good bandwidth efficiency (same as 8-PSK) 
and relatively good noise immunity, the BER of 8-PSK is lower than that of 8-DPSK 
for £= 0°, as expected 

The value of using 8-DPSK versus 8-PSK can be appreciated more fully when 
comparing the performance plots of Figures 4.1 and 4.5. These show that if the 
receiver has a local oscillator phase error restricted to 0 ^ ^ 66.5°, for similar SNR 
values, the BER of the 8-DPSK receiver is lower than that of the 8-PSK receiver. 
Since such phase errors can be typically expected to be small, this is a significant 
advantage. However it must be pointed out that for phase errors of more than 10°, 
both receivers exhibit severely degraded performance that would make both systems 
unsuitable for operation in a practical system. The 8-DPSK receiver does not exhibit 
complete insensitivity to phase errors as encountered in phase measurement receivers. 
Only when a complete phase revesal takes place (i.e., £= 180°) is the performance of 
the 8-DPSK DBDR unaffected, while the 8-PSK receiver produces a nearly useless 
output. 

The angle a set to 22.5° is optimum in minimum BER sense for the 8-PSK 
DBDR. The variations in PE versus a are shown in Figure 4.6, for £= 0° for 8-PSK, 
and a similar plot of PE versus a is plotted for various values of SNR and £= 0° on 
Figure 4.7 for the 8-DSPK case. 
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Figure 4.7 PE Versus a of 8-DPSK DBDR for Various Values of SNR (£= 0°). 
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V. PERFORMANCE ANALYSES OF 16-PSK AND 16-DPSK DBDR'S 
WITH LOCAL OSCILLATOR PHASE ERROR 



In this chapter, DBDRs will be analyzed for 16-PSK and 16-DPSK 
communication systems operating in the presence of a receiver local oscillator phase 
error. The analysis procedure will be similar to that carried out in the previous chapter 
for 8-PSK and 8-DPS K systems. 

The performance of 16-PSK DBDR's has been investigated and the methodology 
for carrying out the performance analyses have been well-documented by Bukofzer 
[Ref 8], In this research, these analyses will be extended to the 16-DPSK system as 
weU as the 16-PSK system modified to operate in the presence of a receiver local 
oscillator phase error. 

A. 16-PSK DBDR PERFORMANCE ANALYSIS IN THE PRESENCE OF 

RECEIVER LOCAL OSCILLATOR PHASE ERROR 

The possible transmitted signals, which are represented as vectors in Figure 1.2.b, 
are mathematically given by Equation 4.1 with i= 1,2,3,.. .,16. The same notation 
introduced in Chapter IV will be used to represent the signals using the orthonormal 
set given by Equation 4.4. Therefore the transmitted signals are expressed as 



S;(t)= (pi(t) cos0j(t) -I- (p^it) sina(t) 
where i= 1,2,... ,16, O^t^T^ and 



(5.2) 



6i(t) 



^ a + (i-l) ;t/8 
p + (i-2) 7t/8 
= < (i+ 1) 7t/8 -P 



i rt/8 -a 



i= 1,5,9,13 
i= 2,6,10,14 
1=3,7,11,15 
1 = 4,8,12,16 



(5.3) 



with a and P as shown on Figure 1.2.b 

Assuming all signals are equally likely to be transmitted, for the 16-PSK case 
Equation 4.7 can be written as 
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16 

Pr{ bit correct } = 1/16 ^ Pr{ bit correct / s.(t) transmitted } (5.4) 

i= 1 

The notation to denote the four bits that make up a data state has been introduced in 
Chapter I. Also the notation introduced by Equation 4.13 is used here, namely 

Pr{ bit correct / Sj(t) } = Pr{ bit correct / Sj(t) transmitted } i= 1,2.., 16 (5.5) 

From Figure 1.2.b, the conditional probabilities for the correct reception of the LB, 
MBL, MBR, RB are given by 



Pr{ LB correct / Sj(t) } = < 



Pr{ r^ •> 0 / S;(t) } i= 1,2,3,4,13,14,15,16 

Pr{r. .<0/s.(t)} 1=5,6,7,8,9,10,11,12 

\ 1 ,1 1 



(5.6) 



f 1.2,7,8,9,10,15,16 

Pr{ MBL correct / Sj(t) } =< (5.7) 

[Pr{ ri/-r2/<0/ Sj(t) } 1=3,4,5,6,11,12,13,14 



Pr{ MBR correct / Sj(t) } = * 



Pr{ RB correct / s^(t) } 



Pr{ 0 / S;(t) } i= 1,4,5,8,9,12,13,16 

Pr{ Vyj^<0 / S;(t) } 1 = 2,3,6,7,10,11,14,15 



|Pr{r2^j>0/Sj(t)} i= 1,2,3, 4,5,6,7,8 

Pr{ r,.<0 / s.(t) } 1=9,10,11,12,13,14,15,16 



(5.8) 



(5.9) 



where the r. . are given by Equation 4.9 and 

K,1 



= 



(5.10) 



The analyses methods which have been given in Chapter IV, can be applied to 
determining the conditional error probabilities for LB, MBL, RB. For MBR, the 
analysis procedure which has been used by Bukofzer [Ref 8], will be utilized. In order 
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to define the conditional probability sets of Equation 5.8, additional steps must be 
completed as demonstrated in the sequel. If Sj(t) transmitted, the MBR is correctly 
recovered if for i= 1,4,5,8,9,12,13,16, 0 and if for i= 2,3,6,7,10,1 1,13,14, < 0. 

Even though the r^ j are conditional Gaussian random variables, because of the highly 
nonlinear form of Equation 5.10, a more productive approach to the determination of 
the conditional p.d.f. of has been introduced (see Reference 8 for details) rather 
than using transformations of Gaussian random variables. By introducing the 
following notation, namely 



j= component of r(t) along <Pj(t) = t cos r\ (5.11) 

r 2 j= component of r(t) along 92^0 ~ C sin q (5.12) 

it can be demonstrated that 

Vih=€'*cos4n (5.13) 

where 

t = length of the vector = V(ri ;)^ "b (f-> :)^ (5.14) 

v\ = vector angle that with respect to the positive x axis (5.15) 

= arctan (r 2 • / rj .) 



Since 1“^ > 0, it is clear that the test involving being greater than or less than zero 
is now equivalent to the test involving cos4q in the same manner. This implies that 
the p.d.f. of T\ must be obtained which, as shown in Reference 8, leads to integrals that 
cannot be evaluated in closed form and therefore must be evaluated via numerical 
integration. 

By utilizing the analysis procedure from Reference 8 and a methodology 
introduced in the previous chapter for the 8-PSK case, the conditional probabilities for 
the correct reception of a bit in the presence of receiver local oscillator phase error is 
given in shorthand form by the following expressions 

Pr{ LB correct / s.(t) } = 1 - Q{ f («,P.2) }, i= 1,2, 3. ..16 (5.16) 

w’here 
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cos (a-c) 


i=l,9 


cos (P-e) 


i = 2,10 


sin (P + e) 


i=3,ll 


sin (a + e) 


i = 4,12 


sin (p-e) 


i=5,13 


sin (a-e) 


i=6,14 


cos (P + e) 


i=7,15 


cos (a + c) 


i=8,16 



and a and p are shown on Figure 1.2.b. Furthermore 



(5.17) 



Pr{ MBL correct / Sj(t) } = Q(-V Es'^^’o Y + sin y) . Q(-VEs/Nq (cos y-sin y) } + 

•+• Q{VEj/Nq (cos y-sin y) }. Q(VEj/Nq (cos y + sin y) } (5.18) 



where 



y = r (a-c) i= 1,5,9,13 
\(p-c) i= 2,6,10,14 



(5.19) 



or 



Pr{ MBL correct / Sj(t) } = 1 - Q{-VEj/Nq (cos y + sin y) . Q{-VEj/Ng (cos y-sin y) } + 

+ Q(VE 3 /Nq (cos y + sin y)} . Q(VEj/Nq (cos y-sin y)} (5.20) 



where 

y=np-£) i=3,7,ll,15 (5.21) 

\ (a-e) i= 4,8,12,16 

For the calculations involving the conditional probability for the MBR, the joint 
p.d.f. of the random variables t| and t has been determined by Bukofzer [Ref. 8] 
namely 

ff (L , H) = |L| f_ _ (L cosH , L sinH) + (5.22) 

i’^2 i/s. 

+ |L1 f*_ ’ * (-L cosH , -L sinH) 

^l,i>*^2,i/Sj 
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which must be integrated to produce the marginal p.d.f. of n- This results in 



(H / Sj) = 1/2 exp(-E 3 /No) + { 1 + V • 



(5.23) 



. { exp[ E 3 /N 0 cos^ (H-£;)] cos(H-^j) . Q[ -^2 E^/Nq (cos(H-£;) ] } 



where 



^; = a(t)-£ i=l,2....16 (5.24) 

where 6 j(t) is given by Equation 5.3 Therefore, the probability that involves 

computing the probability that cos4i]>0. Observe that cos4i]>0 for -7t/8 < rj < 7i;8, 
3JC/8 < n <57t/8, 77t/8< n<97t/8, and 1 l7t/8 < n < 13JI/8. Therefore the calculation of 
the conditional probability of the MBR being correct given that Sj(t) was transmitted 
can be calculated by integrating the conditional p.d.f of r| over the appropriate 
regions, namely 

Pr{ MBR correct / Sj(t) } = Pr{ cos4t]>0 / s^(t) } = 

Jt /8 57t/8 971/8 

- 1 fn/s, + 1 fn/s, 1 'n/s, 

-7t/8 37t/8 77t/8 

137t/8 

+ J f^/ 3 ^ [H/sJ dH 
1171/8 
7t/8 

= f { f^/ 3 ^[H / Sj] + f^^ 3 ^ [(H + 7t/2) / sJ + f^/ 5 ^ [(H + tt) / Sj] + (5.25) 

-7t/8 

+ f^/s^ [(H + 37t/2)/sJ}dH 

By utilizing the appropriate form for the conditional p.d.f. of r| when Sj(t) is assumed 
transmitted, results in 

Pr{ .MBR correct / Sj(t) } = 1/2 exp(-E 3 /Ng) + 
tt/8 

+ V E^/TtNg f { exp[- E 3 /N 0 sin2(H-X)] cos(H-X) . {Q[-V2 E^/Ng (cos(H-X)] - 
-7t/8 
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- Q[J2 Ej,Nq (cos(H + X)] } + { exp[-Ej/NQ cos“(H-X)] sin(H-X) Q[-V2 E^/Ng (sin(H-A)] - 
-Q[vTE7!^(sin(H + X)]} } dH (5.26) 



where 



X = 



a-i-£ 


1=4,8,12,16 




a-£ 


i= 1,5,9,13 


(5.27) 


p-£ 


1=2,6,10,14 




P + £ 


1=3,7,11,15 





Finally, it can be seen that 



Pr{ RB correct / Sj(t) } = Pr{ LB correct / Sj(t) } (5.28) 

so that by using Equations 5.4 and 5.16-5.28 for correct reception probabilities of the 
individual bits, the overall BER can be obtained by using 

Pr( bit in error } = 1 - Pr( bit correct } (5.29) 

This results in the (unconditional) probabilities 

Pr( LB in error } = Pr( RB in error } 

= 1/8 ( Q[ J2 E^/Nq cos(a-c)] + Q[ cos(P-c)] + 

+ Q[ V2 E3/N0 sin(p + c)] + Q[V2E7 Nq sin(a + c)] + 

+ Q[ V2 Ej/Ng sin(P-c)] + Q[V2 E^/Nq cos(P + e)] + 

+ Q [V2 E^/Xq cos(a + e)] + Q [V2 e}Nq sin(a + c)]l (5.30) 



and 

Pr{ .MBL in error } = 1/2 + 1/4 [ Q {-VE^/Nq [cos(a-£) + sin(a-£)]} + 

-H Q {- VEj/Nq [cos(a + £) - sin(a + £)]} + Q {-VEj/Nq [cos(p-£) + sin(p-£)]} + 
-H Q {-VEj/Nq [cos(p-£) - sin(p-£)]} - Q {-VEj/Nq [cos(p + £) + sin(P + £)]} - 
- Q {-VE3/N0 [cos(P + £) - sin(P + £)]} - Q (-VE3/N0 [cos(a + £) + sin(a-H£)]} - 
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- Q {-VEj/Nq [cos(a + £) - sin(a + e)]} ] + 1/2 [ Q {-VE/Nq [cos(a + e) + sin(a + £)]} . 
. Q {-Ve}Nq [cos(a-£) - sin(a-£)]} + Q {-VE^/Nq [cos(P + £) + sin(P + £)]} . 

. Q {-VEj/Ng [cos(p + £) - sin(P + £)]} - Q {VE^/Ng [cos(a-£) + sin(a-£)]} . 

. Q (VEj/Ng [cos(a-£) - sin(a-£)]} - Q {-VEj/Ng [cos(p-£) + sin(p-£)]} . 

. Q {-VeJN^ [cos(P-£) - sin(p-£)]} ] (5.31) 

For the MBR the result appears in integral form. That is 

71/8 

Pr( MBR in error ) = 1 - 1/2 V Ej/TtNg J { exp{-Ej/Ng sin^[H-(a-£)]} cos[H-(a-£)]. 

-7t/8 

. ( 1/2 - Q[-VE7N"g cos(H-(a-£))]} + exp(-Ej/Ng cos^fH-'a-£)]} sin[H.(a-£)] 

. ( 1/2 - Q[-Ve7i^ sin(H-(a-£))]} + expl-E^/Ng sin-[H-(P-£)]} cos[H-(P-£)] 

. { 1/2 - Q[-VE7i^o cos(H-(p-e))]} + expl-E^/Ng cos2[H-(p-£)]) sin[H-(p-£)] 

. { 1/2 - Q[-Ve 7^ sin(H-(P-£))]} - exp{-E^/Ng sin^[H + (p + £)]} cos[H + (P + £)] 

. [ 1/2 - Qi-VE^Ng cos(H + (p +£))]} - exp{-Ej/Ng cos2[H + (P + £)]} sin[H + (p + £)] 

. { 1/2 - Q[-VEj/Ng sin(H + (a + £))]} + exp{-Ej/Ng sin^[H + (a + £)]} cos[H + (a + £)] 
. ( 1/2 - Q[-VEj/Ng cos(H + (a + £)]} + exp{-Ej/Ng cos^[H + (a + £)]} sin[H + (a + £)] 

. ( 1/2 - sin(H + (a + £))]} } dH (5.32) 

The overall BER can now be determined from Equations 5.30,5.31 and 5.32 by 
evaluating the final result for PE, the bit error probability, namely 

PE = 1/4 { Pr{ LB in error } + Pr( MBL in error } 

+ Pr( MBR in error } + Pr{ RB in error } } (5.33) 

It can be seen that Equation 5.33 consists of many terms that by themselves do not 
yield as much insight as the computed values of PE presented in plotted form. 
Therefore plots of BER vs. SN'R for various values of a and £, are given as Figures 5.1 
and 5.2. While a and P can be set independently, from an optimization point of view, 
a + P = 7t/4 has been shown to produce the lowest BER for given values of SNR when 
£ = 0*" [Ref. 8] 
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Figure 5.2 pg Versus SNR for 16-PSK DBDR for Various Values ofe (a= 11.25*'). 
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B. 16 -DPSK DBDR PERFORMANCE ANALYSIS IN THE PRESENCE OF 
RECEIVER LOCAL OSCILLATOR PHASE ERROR 

Similar analyses on DBDR performance have been carried out for 16-DPSK (see 
Chapter IV.B), so that the notation introduced previously will be used here. That is, 

Ck. = Pr { k correct } , Ek = Pr { k in error } (5.34) 

where now k= 1 for LB, k= 2 for MBL, k= 3 for MBR and k= 4 for RB. 

In order to evaluate Pr { LB correct } with differential decoding, we have to 
evaluate 

Pr{ LB correct } = 1/2 [ Pr{ LB correct / LB = 0 } + 

+ { Pr( LB correct / LB = 1 } ] (5.35) 



where 

Pr{ LB correct / LB = 0 ) = (5.36) 

1 1 1 

= Z Z Z Pr( LB correct / LB = 0, MBL = i, MBR = j, RB = k } x 
i=0 j=0 k=0 

X Pr{ MBL = i, MBR = j, RB=k} 

where for i, j, k= 0, 1 

Pr{ MBL = i, MBR = j, RB = k)= 1/8 (5.37) 

In Chapter IV, it has been explained that with DPSK modulation in conduction 
with bit direct detection methods, certain bit errors in the past and presently recovered 
symbol states do not result in errors in the recovered data states. Such bit errors in 
successive symbol states must be considered, because, even though symbol errors may 
occur, the delivered data may still be correct. Therefore, these errors have to be 
accounted for and considered in the computation of the probability of a bit being 
correct. Figure 5.3 show's one example of the situation being discussed. That is, the 
data state (0,0,0,0) is assumed to occur and the symbol state (0,0, 0,0) is assumed to 
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CORRECT RECEPTION COMBINATIONS IN PRESENCE OF ERROR 

INPUT PREVIOUS PRESENT PREVIOUS PRESENT OUTPUT 

TRANSMISSION TRANSMISSION RECEPTION RECEPTION 


0000 0000 0000 0000 


0 0 0 0 


0 0 0 0 


0 0 0 1 


0 0 0 1 


0 0 0 0 


0 0 10 


0 0 10 


0 0 0 0 


0 0 11 


0 0 11 


0 0 0 0 


0 10 0 


0 10 0 


0 0 0 0 


0 10 1 


0 10 1 


0 0 0 0 


0 110 


0 110 


0 0 0 0 


0 111 


0 111 


0 0 0 0 


10 0 0 


10 0 0 


0 0 0 0 


10 0 1 


1 0. 0 1 


0 0 0 0 


10 10 


10 10 


0 0 0 0 


10 11 


10 11 


0 0 0 0 


110 0 


110 0 


0 0 0 0 


110 1 


110 1 


0 0 0 0 


1110 


1110 


0 0 0 0 


1111 


1111 


0 0 0 0 


0 0 0 0 


10 0 0 


0 0 10 


0 0 0 1 


- 1 1 0 


0 0 10 


0 0 10 


110 0 


0 0 10 


0 0 11 


0 111 


0 0 10 


0 10 0 


10 10 


0 0 10 


0 10 1 


0 10 0 


0 0 10 


0 110 


0 0 10 


0 0 10 


0 111 


0 0 0 0 


0 0 10 


10 0 0 


1111 


0 0 10 


10 0 1 


0 0 0 1 


0 0 10 


10 10 


1110 


0 0 10 


10 11 


0 10 1 


0 0 10 


110 0 


110 1 


0 0 10 


110 1 


0 0 11 


0 0 10 


1110 


10 0 1 


0 0 10 


1111 


10 11 


0 0 10 


0 0 0 0 


110 1 


0 10 0 


0 0 0 1 


10 10 


0 10 0 


0 0 10 


10 0 1 


0 10 0 


0 0 11 


0 0 10 


0 10 0 


0 10 0 


1111 


0 10 0 


0 10 1 


10 0 0 


0 10 0 


0 110 


1110 


0 10 0 


0 111 


110 0 


0 10 0 


10 0 0 


0 0 11 


0 10 0 


10 0 1 


0 10 0 


0 10 0 


10 10 


10 11 


0 10 0 


10 11 


0 0 0 0 


0 10 0 


110 0 


0 0 0 1 


0 10 0 


110 1 


0 110 


0 10 0 


1110 


0 10 1 


0 10 0 


1111 


0 111 


0 10 0 


0 0 0 0 


110 0 


0 110 


0 0 0 1 


0 10 0 


0 110 


0 0 10 


1110 


0 110 


0 0 11 


0 110 


0 110 


0 10 0 


10 0 0 


0 110 


0 10 1 


0 0 0 0 


0 110 


0 110 


10 10 


0 110 


0 111 


0 0 10 


0 110 


10 0 0 


110 1 


0 110 


10 0 1 


0 10 1 


0 110 


10 10 


1111 


0 110 


10 11 


0 111 


0 110 


110 0 


10 0 1 


0 110 


110 1 


0 0 0 1 


0 110 



Figure 5.3 Example of Error Patterns for 16-DPSK System Where Data State, 
Previous Transmission, and Present Transmission are 0000 with RB Correctly Recovered. 
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INPUT 


^PREVIOUS PRESENT 

transmission transmission 


PREVIOUS 

RECEPTION 


PRESENT 

RECEPTION 


OUTPUT 






1 


1 


1 


0 


1 


0 


1 


1 


0 


110 






1 


1 


1 


1 


0 


0 


1 


1 


0 


110 






0 


0 


0 


0 


0 


1 


0 


0 


1 


0 0 0 






0 


0 


0 


1 


0 


0 


1 


1 


1 


0 0 0 






0 


0 


1 


0 


0 


0 


0 


0 


1 


0 0 0 






0 


0 


1 


1 


1 


0 


1 


1 


1 


0 0 0 






0 


1 


0 


0 


0 


1 


1 


0 


1 


0 0 0 






0 


1 


0 


1 


0 


0 


0 


1 


1 


0 0 0 






0 


1 


1 


0 


0 


1 


1 


1 


1 


0 0 0 






0 


1 


1 


1 • 


0 


1 


0 


1 


1 


0 0 0 






1 


0 


0 


0 


1 


0 


1 


0 


1 


0 0 0 






1 


0 


0 


1 


1 


1 


0 


1 


1 


0 0 0 






1 


0 


1 


0 


0 


0 


1 


0 


1 


0 0 0 






1 


0 


1 


1 


1 


0 


0 


1 


1 


0 0 0 






1 


1 


0 


0 


1 


0 


0 


0 


I 


0 0 0 






1 


1 


0 


1 


1 


1 


1 


1 


1 


0 0 0 






1 


1 


1 


0 


1 


1 


0 


0 


1 


0 0 0 






1 


1 


L 


L 


1 


L 


1 


0 


1 


0 0 0 






0 


0 


0 


0 


0 


1 


0 


1 


1 


0 1 0 






0 


0 


0 


1 


1 


1 


0 


1 


1 


0 10 






0 


0 


1 


0 


0 


1 


1 


1 


1 


0 10 






0 


0 


1 


1 


1 


1 


1 


1 


1 


0 1 0 






0 


1 


0 


0 


0 


0 


0 


1 


1 


0 10 






0 


1 


0 


1 


1 


0 


0 


1 


1 


0 10 






0 


1 


1 


0 


0 


0 


1 


1 


1 


0 10 






0 


1 


1 


L 


1 


0 


1 


1 


1 


0 10 






1 


0 


0 


0 


0 


1 


0 


0 


1 


0 10 






1 


0 


0 


1 


1 


1 


0 


0 


1 


0 10 






1 


0 


1 


0 


0 


1 


1 


0 


1 


0 10 






1 


0 


1 


1 


1 


1 


1 


0 


1 


0 10 






1 


1 


0 


0 


0 


0 


0 


0 


1 


0 10 






1 


1 


0 


1 


1 


0 


0 


0 


1 


0 10 






1 


1 


1 


0 


0 


0 


1 


0 


1 


0 10 






1 


1 


1 


1 


1 


0 


1 


0 


1 


0 10 






0 


0 


0 


0 


1 


0 


0 


1 


1 


10 0 






0 


0 


0 


L 


1 


0 


0 


0 


1 


10 0 






0 


0 


1 


0 


1 


0 


1 


1 


1 


10 0 






0 


0 


1 


L 


1 


0 


1 


0 


1 


10 0 






0 


1 


0 


0 


1 


1 


0 


1 


1 


10 0 






0 


1 


0 


1 


1 


1 


0 


0 


1 


1 0 0 






0 


1 


1 


0 


1 


1 


1 


1 


1 


10 0 






0 


1 


1 


1 


1 


1 


1 


0 


1 


10 0 






1 


0 


0 


0 
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1 


1 


10 0 
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1 


0 
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0 


1 


10 0 
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0 
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1 


1 


10 0 






1 


0 


1 


1 


0 


0 


1 


0 


1 


10 0 






1 


1 


0 


0 


0 


1 


0 


1 


1 


10 0 






1 


1 


0 


1 


0 


1 


0 


0 


1 


10 0 






1 


1 


1 


0 


0 


1 


1 


1 


1 


10 0 






1 


1 


1 


1 


0 


1 


1 


0 


1 


10 0 






0 


0 


0 


0 


0 


0 


0 


1 


1 


1 1 0 






0 


0 


0 


1 


1 


1 


1 


1 


1 


1 1 0 






0 


0 


1 


0 


0 


1 


0 


1 


1 


110 






0 


0 


1 


1 


1 


1 


1 


0 


1 


110 






0 


1 


0 


0 


0 


0 


1 


1 


1 


110 






0 


1 


0 


1 


1 


1 


0 


1 


1 


110 






0 


1 


1 


0 


1 


0 


1 


1 


1 


110 






0 


1 


1 


1 


1 


0 


0 


1 


1 


110 






1 


0 


0 


0 


0 


1 


1 


0 


1 


1 1 0 






1 


0 


0 


1 


1 


0 


0 


0 


1 


110 






1 
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0 
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1 


110 






1 


0 
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1 
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1 


0 


0 


1 


110 






1 


1 


0 


0 


0 


1 


0 


0 


1 


110 






1 
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1 


1 
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110 
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Figure 5.3 Continued. 
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have been previously transmitted. This means that according to the encoder matrix 
shown in Chapter II, the symbol state (0,0, 0,0) must now be generated for 
transmission. Figure 5.3 shows the receptions of the symbol states which produce 
recovered data states in which the RB is correctly delivered. Observe that there are 128 
such combinations which yield 0 for the RB. Therefore notation has to be introduced 
to account for error patterns of the above type, namely Pr( RB correct, (g,h,p,q) / 

LB = m,MBL = i,MBR = j,RB = k }, where g,h,p,q= 0,1. The meaning of this 
expression is the probability of the RB being correctly recovered when the data state is 
(m,j,i,k) and the delivered data state is (g,h,p,q) where q = 0 in this case (otherwise the 
RB could not be correct due to the transmission of the RB = 0). Therefore 

Pr( RB correct / LB = 0, MBL = 0 ^'^BR = 0, RB = 0 } = 

1 1 1 

= Z Z Z RB correct, (g,h,p,0) / LB = 0, MBL = 0, MBR = 0, RB = 0 } (5.38) 

g = 0 h = 0 p = 0 

Since each element of the symbol state is independently recovered, we see from Figure 
5.3 that 

Pr{ RB correct,(0,0,0,0) / LB = 0,MBL = 0,MBR = 0,RB = 0 } = Cl^ C2^ C3^ C4^ + 

+ Cl^ C2^ E3^ C4^+ Cl^ C2^ E3^ E4^+ Cl^ E2^ C3^ E4^+ Cl^ E2^ E3^ C4^ + 

+ Cl- E2^ E3^ E4^+ El^ + C2^ C3^ C4^+ El^ C2^ C3^ E4^ El^ C2^ E3- C4^ + 

+ El^ C2^ E3^ E4^+ El^ E2^ C3^ C4^+ El^ E2^ E3^ C4^+ El^ E2^ E3^ E4^ + 

+ Cl- C2^ C3^ E4^ (5.39) 

However, from Equations 5.34, and 5.30, it is shown that Cl= C4 and El= E4. 
Therefore, 

Pr( RB correct, (0,0,0,0) / LB = 0, MBL = 0, MBR=0, RB = 0 } = 

= Cl^ C2^ C3-+ 2 Cl^ C2^ C3^ El^ + C2^ C3^ El^ E2‘ 

+ 2 Cl-^ C2^ E3^ + 2 Cl^ C2^ El^ E3^+ Cl^ C3^ E2^ + El^ E2^ E3^ 

2 Cl- C3^ El^ E2^ + C2^ El^ E2^ E3^ + Cl^ El^ E2^ E3^ + 

^ C3^ El^ E2^ + Cl^ C3^ El'^ (5.40) 
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A similar expression must be obtained for other combinations of values of g,h and p as 
well as the condition that RB is 1. This produces a total of 256 possibilities of received 
symbol states that result in a correct RB. Computer generated outputs were used to 
identify the equivalent error patterns for other cases involving LB=1, MBL=1 and 
MBR= 1. These results are shown in Appendix C. Similar analyses must be carried 
out in order to determine the unconditional probabilities Pr{LB correct}, Pr{MBL 
correct) as well as Pr{MBR correct). Such analyses demonstrate due to inherent 
symmetries that 



Pr{LB correct) = Pr{MBL correct) = Pr{MBR correct) = Pr{RB correct) = PE 



so that using Equation 5.38 along with similar forms obtained for k~ I and 1= 1 yield 



PE = 1 - 1/2 { 2 Cl-^ C2^ C3^ + Cl^ C2^ C3 E3 + 3 Cl^ C2^ C3 E2 + 

+ 2 Cl^ C2 C3 E2 E3 H- 6 Cl^ C2^ C3^ El + 4 Cl^ C2^ C3 El E3 + 

+ 12 Cl^ C2 C3- C3 El E2 + 8 Cl^ C2 C3 C3 El E2 + 6 Cl^ C2^ El E3^ + 

+ 6 Cl^ C2- C3" El + Cl** C3 E2" E3 + 4 Cl^ C3 El E2^ E3 + 

+ 2 C1-* El^ E3^ + Cl^ C2^ E3^ + 3 Cl^ C2 E2 E3^ + 2 C2^ C3^ El‘^ + 

+ Cl'^ C3^ E2^ + 10 Cl^ C2^ C3^ El^ + 10 Cl^ C2^ El^ E3' + 

+ 10 Cl^ C3^ El^ E2^ + 8 Cl^ El^ E2^ E3^ + 18 Cl^ C2 C3^ El^ E2^ + 

+ 18 Cl- C2 El^ E2 E3^ + 6 Cl^ El E2^ E3^ + 12 Cl^ C2 El E2 E3^ + 

+ 12 Cl^ C2 C3 El^ E2 E3 + 6 Cl^ C2^ C3 El^ E3 + 5 Cl^ C3 El^ E2^ E3 + 

+ 6 Cl C2^ C3^ El^ + 6 Cl C3^ El^ E2^ + 6 Cl C2^ El^ E2^+ C2^ El*^ E3^ 
+ 12 Cl C2 C3^ El^ E2 + 12 Cl C2 El^ E2 E3^ + 6 Cl El^ E2^ E3^ + 

+ 8 Cl C2 C3 El^ E2 E3 + 4 Cl C3 El^ E2^ E3 + 4 Cl C2^ C3 El^ E3 + 

+ C3^ El^ E2^+ 2 El^ E2^ E3^ + C2^ C3 El'^ E3 + C3 El'^ E2^ E3 + 

H- 3 C2 El-^ E2 E3^ + 2 C2 C3 El"^ E2 E3 + 2 C2 C3^ El'* E2 ) (5.42) 



In Figure 5.4, PE (given by Equation 5.42) is plotted for various values of a. Also 
Figure 5.5 includes BER for various values of e when a= 11.25° 



(5.41) 
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C. NUMERICAL RESULTS AND DISCUSSION OF RESULTS FOR 16-PSK 
AND 16-DPSK DBDR'S WITH RECEIVER LOCAL OSCILLATOR PHASE 
ERROR 

In this chapter 16-PSK and 16-DPSK DBDRs have been analyzed when 
operating in the presence of a receiver local oscillator phase error. The performance 
results which have been obtained for the 16-ary case, are similar to those obtain for 
the 8-ary case. Hence, instead of repeating the conclusions similar to these presented 
for the 8-ary case, conclusions that are unique to the 16-ary case will be presented in 
this part. Main differences between 8-ary and 16-ary cases which are observed via 
analyses and numerical evaluations are 

1. The optimum angles in minimum BER sense, a and P between the symbols 
becomes 11.25° and 33.75° respectively for 16-PSK case, as well as for 
1 6-DSPK case, with a + p = 45° 

2. The usage of 16-PSK provides better performance than 16-DPSK for 8^33.5°, 
however the actual BER values are too high for use in a practical system when 
8 greater than a few degrees. 

3. The 16-DPSK system performs better than the 16-PSK system for 8>33.5°. 
However again the actual values of PE under these circumstances is too high 
for a practical system. 

4. The 16-DPSK system is insensitive to phase errors of 7C radians. 

Performance comparisons between the 16-PSK and 16-DPSK DBDRs are plotted 

on Figure 5.6, and BERs are plotted on Figures 5.7 and 5.8 respectively for various 
values of a with 8=0° 
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SNR(DB) 



Figure 5.6 PE Versus SNR Comparison Plot for 16-PSK and 16-DPSK DBDR 

(«= 11.25^E= 0«). 
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Figure 5.7 PE Versus a for 16-PSK DBDR for Various Values of SNR (€= 0°). 
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Figure 5.8 PE Versus a for 16-DPSK DBDR for Various Values of SiS'R (e= O'^), 
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VI, CONCLUSIONS 



In this research the use of DBDRs had been analyzed for M-PSK and MDPSK 
communication systems. The main assumptions in the analyses are a Gray coding 
mapping for the symbols and an Additive White Gaussian Noise channel model. It is 
funher assumed that the receiver has a known phase error in its local oscillator and its 
effect on receiver performance is analyzed for the the special cases in which M = 8 and 
M=16. 

In addition to the analytical results, graphical results are presented in the form of 
plots of variations in PE for various values of phase angle a, SNR, and local oscillator 
phase error £. 

Optimum conventional M-PSK and M-DPSK receivers are designed to detect 
symbols while DBDRs are designed to detect bits. However, local oscillator phase 
errors degrade the performance of .M-PSK and M-DPSK conventional receivers. Part 
of the work undertaken was intended to investigate the effect of such phase errors on 
DBDR's. The variations of PE versus £ are plotted in Figure 6.1 for 8-PSK and 
8-DPSK, and in Figure 6.2 for 16-PSK and 16-DPSK. If £ (receiver local oscillator 
phase error) can be constrained to a small interval, namely 0°^£^30°, the xM-PSK 
system provides better performance than the M-DPSK system. For values of £ outside 
this range, both systems have high error rates that make the receivers unsuitable for 
practical applications. However, for complete phase reversal, i.e., £= 180°, the M- 
DPSK system suffers no performance penalty whereas the M-PSK system can not 
properly operate. 

The performance degradation due to a phase error has been analyzed also from 
a different point of view. That is, when £ ^ 0°, there is an increase in PE for both the 
M-PSK and M-DPSK receivers. The ratio of the phase error is given by 

PER=PEg/PE^ (6.1) 

where 

PEg = Probability of bit error @ £ = 0° 

PEj = Probability of bit error @ £ * 0° 
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PER has been analyzed under various conditions for M-PSK and M-DPSK. The 
results have been plotted in Figures 6.3and 6.5 for M = 8, and in Figures 6. 6-6. 8 for 
M=16, as a function of SXR for different values of s. The plots demonstrate 
significant differences in the performance degradation of the M = 8 versus M= 16 case. 
However, for high SNR values, these differences tend to vanish. 

The analysis results lead to the following conclusions 

1. By using DBDRs, the error probability of each bit in the symbol state can be 
calculated independently and more importantly these error probabilities can be 
varied by changing the angle between symbols. 

2. For optimum BER, a should be set to 22.5*^ for 8-PSK and 8-DPSK systems, 
while 16-PSK and 1 6-DPS K systems require a= 11.25® and P + a = 45®. 

3. The solution for the phase ambiguity problem cannot be achieved by using M- 
DPSK DBDR's, unless the phase ambiguity is 180®. Therefore, if the perfect 
coherency conditions cannot be achieved, the M-PSK system should be used 
whenever DBDR's are used. 

4. When c= 0®, the DPSK system exhibits good performance. However, if the 
£= 0® condition can be guaranteed there is no need to use the differential 
encoding scheme. 
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Figure 6.3 PER (£ = 0°/e= 5°) Versus SNR for 8-PSK 
andS-DPSK DBDR's (a= 22.5'^). 




Figure 6.4 PER (£ = 0®;c= 22.5'") Versus SNR for 8-PSK 
andS-DPSK DBDR's (a= 22.5'"). 
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Figure 6.5 PER (s = 0°/c = 45°) Versus SNR for 8-PSK 
andS-DPSK DBDR's (a= 22.5°). 




Figure 6.6 PER (€ = 0°'€=5°) Versus SNR for I6-PSK 
and 16-DPSK DBDR's (a = 11.25°). 
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Figure 6.7 PER (£ = 0^^, £= 11.25*^) Versus SNR for 16-PSK 
and 16-DPSK DBDR's (a= 11.2i°). 




Figure 6.8 PER (£ = 0®'£ = 45^) Versus SNR for 16-PSK 
and 16-DPSK DBDR's (a= 11. 25'^). 
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APPENDIX A 

DBDR INPUT/OUTPUT GENERATION PROGRAM (BELIZ) 



1. USER GUIDE 

This program generates all the possible input / output combinations in an M- 
DPSK communication system. The program runs interactively and sends outputs to 
the files 

1. FILE FT08F001 

2. FILE RCDATA 

3. FILE XMDATA 

4. FILE R 

5. FILE ERR. 

These files are used for information transfer between subroutines but are not required 
for program execution. It is possible to generate input / output combinations for 
8-DPSK and 16-DPSK by using this program, while for 32-DPSK or 64-DPSK 
systems, dimensions and format statements must be changed accordingly, and a larger 
character variable set in the subroutine XMTR must be created. The program accepts 
only the integer equivalents of your data or symbol states. An example run of the 
program is shown below 

run beliz 

EXECUTION BEGINS... 

PLEASE SELECT ONE OF THEM 

(1) IF YOU HAVE TRANSMITTER 

(2) IF YOU HAVE RECEIVER 



1 

PLEASE SELECT ONE OF THEM 

(1) IF YOUR SYSTEM 8-DPSK 

(2) IF YOUR SYSTEM 16-DPSK 

1 

ENTER THE INTEGER EQUIVALENTS OF SYMBOL STATES 
IN THE ORDER OF COUNTER CLOCK WISE 
ENTER THE STATE # 1 



7 

ENTER THE STATE # 2 



5 

ENTER THE STATE # 3 



1 

ENTER THE STATE # 4 
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3 



5 



ENTER THE STATE # 



2 

ENTER THE STATE # 6 



0 

ENTER THE STATE # 7 



4 

ENTER THE STATE # 8 



6 

ENTER THE INTEGER EQUIVALENTS OF THE 
PHASE ANGLES OF YOUR CODING SYSTEM 
1*ALPHA= ? 



0 

1 

3 

2 

6 

7 

5 



2*ALPHA= ? 

3*ALPHA= ? 
4*ALPHA= ? 

5*ALPHA= ? 
6*APLHA= ? 
7*ALPHA= ? 
8*ALPHA= ? 



4 

PLEASE ENTER THE SELECTION CORRECTLY 

ENTER 1 IF YOU NEED BINARY FORM OF SYMBOL AND 
DATA STATES INCLUDING SYMBOL OUTPUTS 

ENTER 2 IF YOU NEED DECIMAL FORM OF SYMBOL AND 
DATA STATES INCLUDING SYMBOL OUTPUTS 

ENTER 3 IF YOU WANT TO TURN BACK MAIN PROGRAM. 



1 

PLEASE ENTER THE SELECTION CORRECTLY 

ENTER 1 IF YOU NEED BINARY FORM OF SYMBOL AND 
DATA STATES INCLUDING SYMBOL OUTPUTS 

ENTER 2 IF YOU NEED DECIMAL FORM OF SYMBOL AND 
DATA STATES INCLUDING SYMBOL OUTPUTS 

ENTER 3 IF YOU WANT TO TURN BACK MAIN PROGRAM. 



3 



ANALYSES RESULTS WAS SENT TO THE FILE FT08F001 
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1 



SOURCE CODE 



THIS PROGRAM SIMULATES A M-DPSK COMMUNICATION SYSTEM 
PROGRAM RUNS TOTATLY INTERACTIVE THEN SENDS TO OUTPUTS FILES; 
FILE FT08F001,FILE RCDATA,FILE XMDATA,FILE R,FILE ERR. 

THESE FILES ARE USED FOR INFORMATION TRANSFER BETWEEN 
SUBROUTINES BUT ARE NOT REQUIRED FOR PROGRAM EXECUTION. 

YOU CAN SIMULATE 8-DPSK AND 16-DPSK BY USING THIS PROGRAM 
FOR 32-DPSK OR 64-DPSK SYSTEMS YOU HAVE TO CHANGE DIMENSIONS 
AND FORMAT STATEMENTS ACCORDINGLY , ALSO YOU HAVE TO EXPEND 
CHARACTER VARIABLE SET IN THE SUBROUTINE XMTR. 



OPERATION SPECIFICATIONS. 



THE SYMBOL STATES REQUESTEd DEPEND ON YOUR SYSTEM AND CAN 
BE ANY SEQUENCE. MAIN THING IS; YOU HAVE TO SELECT YOUR REFERENCE 
THEN ENTER THEM BY ONE STEP INCREMENTS , ALSO SAME SPECIFICATION 
EXISTS FOR YOUR ANGLES . PROGRAM ONLY EXCEPTS THE INTEGER EQUIVALENTS 
OF YOUR STATE OR ANGLE CODES. 

:cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 



VARIABLE DECLERATION 

INTEGER N,L, ORDER, TM(64, 64), PS(64) ,TRDATA(256) ,M,KK,K,THETA(64) 
INTEGER MAT^4) ,KL,RMAT(64,64) ,TR(64) ,TH(64,64) ,RCVMAT(64,64) ,NN 
INTEGER QQU,QU,SQU,LD,IS,NT(1,8),KM(1,8),KN(1,8) 

INTEGER RCV(64,64) ,T(256) ,TD(256) ,NEXTSTU56) ,PRESTA(256) ,FL 
INTEGER A,B,Y,KD,KYU,8) ,KH(1,8) ,KJ(1,8) ,KBU »8) 

CHARACTER=^1 DATA (64, 64) 

CHARACTER*16 SET 



*** CALLING SUBROUITINES *** 

CALL INFO (ORDER, PS, THETA, NN,N,FL) 

IF (FL.EQ.l)THEN 
CALL XMTR ( ORDER, FL) 

CALL XMSTA (THETA,TM, PS , ORDER) 

CALL XMDAT ( TM , ORDER , THETA , NEXTST , PS , T , TD ) 

CALL CONVER (N, NEXTST, NN, ORDER, FL) 

WRITE(6,1) 

1 FORMAT(T9,' ANALYSES RESULTS WAS SENT TO THE FILE FT08F001') 
ELSE 

CALL RCVR (ORDER) 

CALL RCVSTA (PS , ORDER, THETA, RCVMAT) 

CALL RCVDAT ( TR , PS , ORDER , NEXTST , RCVMAT ) 

CALL CONVER (N, NEXTST, NN, ORDER, FL) 

CALL ERROR ( FL , THETA , RCVMAT , PS , ORDER , TR , TM , T , NEXTST , TD ) 
WRITE(6,2) 

2 FORMAT(T9,' ANALYSES RESULTS WAS SENT TO THE FILE FT08F001',/, 
*T7, 'ERROR ANALYSES RESULTS WAS SENT TO THE FILE ERR') 

END IF 

STOP 

END 



SUBROUTINE INFO (ORDER , PS , THETA , NN , N , FL ) 

THIS SUBROUTINE PROVIDES INFORMATIONS ABOUT OUR SYSTEM 

INTEGER N.L, ORDER, TM(64, 64) ,PS(64) ,THETA(64) ,FL 
CHARACTER=^1 DATA (64, 64) 

WRITE (6,*) ' PLEASE SELECT ONE OF THEM' 

WRITE (6,*) '(1) IF YOU HAVE TRANSMITTER' 

WRITE (6 ,*) '(2) IF YOU HAVE RECEIVER' 

READ (6,’*') FL 

WRITE (6,^) ' PLEASE SELECT ONE OF THEM' 

WRITE (6,*) '(1) IF YOUR SYSTEM 8-DPSK' 
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WRITE (6,*) '(2) IF YOUR SYSTEM 16-DPSK' 

READ {6,h L 

IF(L.EO.l) THEN 

ORDER=S 

NN=63 

N=6 

END IF 

IF(L.EQ.2) THEN 

0RDER=16 

N=a 

NN=255 

END IF 

WRITE (6,*) 'ENTER THE INTEGER EQUIVALENTS OF SYMBOL STATES' 
WRITE (6,*) 'IN THE ORDER OF COUNTER CLOCK WISE' 

DO 11 K =1, ORDER 
WRITE (6, 20) K 

20 FORMAT U OX, 'ENTER THE STATE #',I3) 

READ(6,*) PS(K) 

11 CONTINUE 

WRITE (6,*) 'ENTER THE INTEGER EQUIVALENTS OF THE ' 

WRITE (6,*) 'PHASE ANGLES OF YOUR CODING SYSTEM ' 

DO 33 L=l, ORDER 
WRITE (6,10) L 

10 FORMAT aOX,I3, '*ALPHA= ?') 

READ (6,*) THETA (L) 

33 CONTINUE 
RETURN 
END 
C 

SUBROUTINE XMTR ( ORDER, FL) 



THIS SUBROUTINE PROVIDES GENERAL TRANSMITER MATRIX FOR 
M-DPSK SYSTEM. 

ITS OUTPUT CAN BE SEEN ON EITHER FILE FT08F001 OR FILE 
XMDATA. 

INTEGER N.L, ORDER, FL 
CHARACTER^ie SET 
SET(1:16)= ' ABCDEFGHIJKLMNOP ' 

OPEN ( UNI T=11,FILE=' XMDATA' , STATUS= ' OLD ' ) 

IF (FL.EQ.l) THEN 
WRITE (8, 81 ) ORDER 

81 FORMAT(///,T15, 'TRANSMITTER'S ENCODER MATRIX ANALYSIS 

* FOR' ,13, '-DPSK' 

*,/,60('-'),//,T5, 'GENERAL TRANSMITTER MATRIX ',/ ,T2 , 31 ('-') ) 
ELSE 

WRITE (8,88) 

88 FORMAT (///,T1 5, 'TRANSMITTER INFORMATION FOR THE GIVEN RECEIVER' 

*,/,72( '-' ),//,T5, 'GENERAL TRANSMITTER MATRIX' ,/ ,T2,31( ' - ' ) ) 

END IF 

DO 3 N=l, ORDER 
IF (N.EQ.l) THEN 
WRITE (8,*) SET(N:ORDER) 

WRITE (11,^) SET (N: ORDER) 

ELSE 

WRITE (8,*) SET (N -.ORDER) .SET (1 :N-1) 

WRITE (11,^) SET(N:ORDER),SET(l;N-l) 

END IF 

3 CONTINUE 

CLOSE (11) 

RETURN 

END 

SUBROUTINE XMSTA( THETA , TM , PS , ORDER) 

THIS SUBROUTINE CONVERTS GENERAL XMTR MATRIX IN TO OUR SPECIAL 
XMTR SYMBOL MATRIX. 

GENERAL XMTR MATRIX IS READ FROM FILE XMDATA. 

INTEGER TM ( 64 , 64 ) , PS ( 64 ) , THETA ( 64 ) , ORDER 
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noon o nonnn 



99 



10 



43 



20 



6 

4 



83 

11 

77 

93 

12 

3 



C 



CHARACTER*! DATA (64, 64) 

0PEN(UNIT=11,FILE=‘XMDATA' , STATUS ='0LD' ) 

IF ( ORDER. EQ. 8) THEN 

WRITE(8, 99) (THETA(I) ,1=1, ORDER) 

FORMAT (///,T5, 'TRANSMITTER MATRIX ',/, 31 ('-'),// , 

*' SYMBOL STATE (COLUMNS) / DATA STATE (ROWS)' 
*,//,4X,3(I3h/35('-')) 

READ (11,10) ((DATA(L,N),N=1,8),L=1,8) 

FORMAT (1X,8A1) 

ELSE 

WRITE (8,43) (THETA ( I ) , 1=1 , ORDER) 

FORMAT (///,T5, 'TRANSMITTER MATRIX ',/, 31 ('-'),// , 

*' SYMBOL STATE (COLUMNS) / DATA STATE (ROWS) 
*,//,4X,16(I3),/ ,70('-'h 
READ (11,20) ((DATA(L,N),N=1,16),L=1,16) 

FORMAT (1X,16A1) 

END IF 

DO 3 K=l, ORDER 
DO 4 N=l, ORDER 
DO 6 M= 0, ORDER 

IF(ICHAR(DATA(K,N)) .EQ.ICHAR(DATA(l,ORDER-M)))THEN 
TM(K,N)=PS(ORDER-M) 

END IF 
CONTINUE 
CONTINUE 

IF ( ORDER. EQ. 8) THEN 
WRITE(11,11) (TM(K,N),N=1, ORDER) 

WRITE(8,83) PS(K), (TM(K,N) ,N=1, ORDER) 

FORMAT(I3, ' ' ,8(13)) 

FORI-IAT (812) 

ELSE 

WRITE(11,12) (TM(K,N) ,N=1, ORDER) 

WRITE (8, 93) PS(K), (TM(K,N) ,N=1, ORDER) 

F0RMAT(I3, ' 1 ' ,16(13)) 

FORMAT (1613) 

END IF 
CONTINUE 
RETURN 
END 

SUBROUT INE XMDAT ( TM , ORDER , THETA , NEXTST , PS , T , TD ) 



THIS SUBROUTINE FINDS THE NEXT SYMBOL STATE/S BY DETERMINING 
XMTR MATRIX CONSTRUCTED BY SUB. XMSTA 



INTEGER TM(64,64) ,PS(64) ,TRDATA(256) ,KK,K,THETA(64) ,NN 
INTEGER ORDER , T ( 256 ) , TD ( 2 56 ) , NEXTST (256), PRESTA ( 256 ) 

DO 5 1=1, ORDER 
T(PS(I))=I 
TD(THETA(I))=I 
5 CONTINUE 
KK=0 

DO 7 K=0, ORDER-1 
DO 8 M=0, ORDER- 1 
KK=KK+1 

NEXTST(KK)=TM(T(K) ,TD(M)) 

8 CONTINUE 
7 CONTINUE 
RETURN 
END 

SUBROUTINE RCVR (ORDER) 

THIS SUBROUTINE PREPARES GENERAL RECEIVER MATRIX AND 
SENDS OUTPUTS TO THE "FILE RCDATA" . 

INTEGER ORDER , MAT ( 64 ) , KL , RMAT (64,64) 

OPEN(UNIT=2,FILE=' RCDATA' , STATUS= ' OLD ' ) 

WRITE (8,81)ORDER 
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oooooo 



81 



3 



44 



66 



77 

55 



10 

15 



11 

16 

88 



C 



FORMAT (///,T1 5, 'RECEIVER'S DECODER MATRIX ANALYSIS 
* FOR' ,13, ' -DPSK' ,/, 

'-') ,//,T5, 'GENERAL RECEIVER MATRIX' ,/ ,T2 , 31 (' - 
DO 3 1=0, ORDER- 1 

mat(ih-iHi 



CONTINUE 
DO 44 K=l, ORDER 
RMAT(1,K)=MAT(K) 

CONTINUE 
DO 55 L=2, ORDER 
DO 66 M=1,L-1 

RMAT(L,M)=MAT(ORDER-L+M+l) 
CONTINUE 
DO 77 N=L, ORDER 
RMAT ( L , N ) =MAT ( N-L+ 1 ) 
CONTINUE 
CONTINUE 

DO 88 NK=1, ORDER 
IF ( ORDER. EQ. 8) THEN 



WRITE (2,10 
WRITE (8,15 
FORMAT (8(12) 



(RMAT(NK,KL 

;rmat(nk,kl 



),KL=1, 

Kkl=i, 



=1, order) 

ORDER) 



F0RMAT(8X,8(I2)) 
ELSE 
WRITE 
WRITE 



( 2 , 1 1 ) ( RMAT ( NK , KL ) , KL= 1 , ORDER ) 
(8,16) (RNAT(NK,KL) ,KL=1, ORDER) 
FORMAT(16(I2)) 

F0RMAT(8X, 16(12)) 

END IF 
CONTINUE 
CLOSE(2) 

RETURN 

END 



)) 



SUBROUTINE RCVSTA ( PS , ORDER , THETA , RCVMAT ) 



THIS SUBROUTINE PREPARES SPECIFIC RECEIVER MATRIX UNDER 
DEFINED STATES.it READS GENERAL RECEIVER MATRIX FORM FILE 
RCDATA WHICH HAD BEEN CREATED BY THE SUBROUTINE RCVDAT. 
OUTPUT ARE SENT TO THE FILE R. 



INTEGER TH(64,64) , ORDER, PS ( 64 ), THETA (64) , RCVMAT (64, 64) 
OPEN (UNIT=2,FILE=' RCDATA' , STATUS= ' OLD ' ) 
OPEN(UNIT=20,FILE='R' ,STATUS= ' OLD ' ) 

IF ( ORDER. EQ. 8) THEN 

WRITE ( 8 , 99) (PS (L ) , L=1 , ORDER) , (THETA( I ) , 1=1 , ORDER) 

99 FORMAT(///,T5, 'RECEIVER MATRIX IS ',/, 31 ('-'),// , 

^'SYMBOL STATE (COLUMNS) / DATA STATE (ROWS) ',// ,4X, 8(13) ) 
ELSE 

WRITE(8, 199) (PS(L) ,L=1 , ORDER) 

199 FORMAT (///,T5, 'RECEIVER MATRIX IS ',/, 31 ('-'),// , 

*' SYMBOL STATE (COLUMNS) / DATA STATE (ROWS)' 

*,//,4X, 16(13)) 

END IF 

DO 1 1=1, ORDER 
IF ( ORDER. EQ. 8) THEN 
READ(2,11) (TH(I,N) ,N=1, ORDER) 

11 FORMAT (81 2) 

ELSE 

READ(2,12) (TH(I,N),N=1, ORDER) 

12 FORMAT (1612) 

END IF 

DO 2 K=l, ORDER 
RCVMAT(I,K)=THETA(TH(I,K)+1) 

2 CONTINUE 

IF ( ORDER. EQ. 8) THEN 

WRITE (20,13) (RCVMAT ( I , N ) , N= 1 , ORDER ) 

13 FORMAT (81 2) 

WRITE (8, 83) PS(I) , (RCVMAT ( I, N) ,N=1, ORDER) 

83 F0RMAT(I3, ' 1 ' ,8(13)) 
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ELSE 

WRITERS, 93) PS(I) , (RCVMAT(I,N) ,N=1, ORDER) 

93 FORMAT 13, ■ 1 M6(I3)) 

WRITE(20,14) (RCVMAT(I,N) ,N=1, ORDER) 

14 FORMAT (161 2) 

END IF 

1 CONTINUE 
CLOSE(20) 

RETURN 

END 

C 

SUBROUTINE RCVDAT ( TR , PS , ORDER , NEXTST , RCVMAT ) 

THIS SUBROUTINE READS SPECIFIC RECEIVER MATRIX WHICH HAD 
BEEN CREATED BY SUBROUTINE RCVSTA,THEN SEARCHS FOR THE RCVR 
OUTPUT. 



INTEGER RCVMAT ( 64 , 64 ) , PS ( 64 ) , KK , ORDER 
INTEGER TR(64) ,NEXTST(256) 
OPEN(UNIT=20,FILE='R' ,STATUS= ' OLD ' ) 

DO 23 LL=1, ORDER 
IF (0RDER.E0.3)THEN 

READ (20,24)(RCVMAT(LL,KM) ,KM=1, ORDER) 
24 FORMAT(8(I2)) 

ELSE 

READ (20,94) (RCVMAT (LL, KM) ,KM=1, ORDER) 
94 F0RI-IAT(16(I2)) 

END IF 

23 CONTINUE 

DO 2 1=1, ORDER 
M=PS ( I ) 

TR(M)=I 

2 CONTINUE 

KK=0 

DO 3 M=0, ORDER- 1 
DO 4 K=0, ORDER- 1 
KK=KK+1 

NEXTST (KK)=RCVMAT(TR(M) ,TR(K)) 

4 CONTINUE 

3 CONTINUE 

RETURN 

END 



SUBROUTINE CONVER ( N , NEXTST , NN , ORDER , FL ) 



INTEGER N,S,QU,SU,LLL,LLK(50) ,MMM,NN,FL,MM,M,K(9,1024) 
INTEGER MW,NT(256,4) , KZ, ORDER, NEXTST(256) , ST, QQU,KD 



PLEASE ENTER THE SELECTION CORRECTLY' 

' ENTER 1 IF YOU NEED BINARY FORM OF SYMBOL AND 
'DATA STATES INCLUDING SYMBOL OUTPUTS' 

' ENTER 2 IF YOU NEED DECIMAL FORM OF SYMBOL AND 
'DATA STATES INCLUDING SYMBOL OUTPUTS' 

' ENTER 3 IF YOU WANT TO TURN BACK MAIN PROGRAM. ' 



300 WRITE (6,*) 

WRITE (6,- 
WRITE 
WRITE 
WRITE 
WRITE 
READ ( 

IF (S.EQ.l 
IF (S.EQ.2 
IF (S.EQ.3 
IF (S.NE.l 
200 IF (FL.EQ.l) 

WRITE (8, 6) 

6 FORMAT (//,8X, 'DATA STATE' ,5X, 'SYMBOL STATE' ,5X, 'SYMBOL' ) 
ELSE 

WRITE (8, 21) 

21 FORMAT (//,8X,'DATA STATE ', 5X, ' SYMBOL STATE ', 5X, ' SYMBOL ' ) 
END IF 
KK=0 



GOTO 100 
GOTO 200 
GOTO 400 

.0R.S.NE.2.0R.S.NE.3) GOTO 300 
THEN 



93 



noon 



9 

8 

7 

100 

47 

77 



33 

44 



22 



45 



55 



46 

23 

400 



DO 7 KZ=0, ORDER- 1 
DO 8 M=0, ORDER- 1 
KK=KK+1 

WRITE(8,9) KZ,M, NEXTST(KK) 

FORMAT (10X,I3,15X,I3,15X,I3) 

CONTINUE 
CONTINUE 
GOTO 300 

IF (FL.EQ.l) THEN 
WRITE (8, 47) 

FORMAT (//,7X,‘DATA STATE ',1 OX, ' SYMBOL STATE ',10X, 

*' SYMBOL' ,/,72('-')) 

ELSE 

WRITE (8, 77) 

FORMAT (//,7X, 'SYMBOL STATE ', lOX, ' DATA STATE ',10X, 
^'SYMBOL' ,/,72('-')) 

ENDIF 

IF(ORDER.EQ.8)THEN 

KD=3 

ELSE 

KD=4 

ENDIF 

DO 44 L=1,NN+1 
QOU=NEXTST(L) 

D5 33 SR=KD,1,-1 
NT(L,SR)=MOD(QQU,2) 

QQU=QQU/2 
CONTINUE 
CONTINUE 
DO 23 MW=0,NN 
QU=MW 

DO 22 LLL=N,1,-1 
K(MW,LLL)=MOD(QU,2) 

QU=QU/2 

CONTINUE 

IF ( ORDER. EQ. 8) THEN 

WRITE(8,45) MW,(K(MW,LK) ,LK=1,N) , (NT(MW+1 ,ST) ,ST=1 ,KD) 
FORMAT (I3,2X,3(2X,I2,2X) ,3X,3(2X,I2,2X),3X,3UX,I2,2X)) 
ELSE 

WRITE (8, 55) MW, (K(MW,LK) ,LK=1,N) , (NT(MW+1 ,ST) , ST=1 ,KD) 
FORMAT (I3,2X,4(2X,I2,1X) , 3X,4(2X, 12 , IX) , 3X,4 (2X, 12 , IX) ) 



ENDIF 
WRITE 
FORMAT 
CONTINUE 
GOTO 300 
RETURN 
END 



8,46) 

72('-')) 



SUBROUTINE ERROR ( FL , THETA , RCVMAT , PS , ORDER , TR , TM , T , NEXTST , TD ) 
THIS SUBROUTINE LOCATES UNDETECTED ERROR COMBINATIONS. 



INTEGER QQU,QU, SOU, PS (64), ORDER, RCVMAT (64, 64) ,N,TR(64) ,TM(64,64) 
INTEGER TD(256) ,T(256),A,B,Y,KD,KN(1,8),KJ(1,8) ,KB(l,a) ,KY(1,8) 
INTEGER FL,IS,NTU,8) ,KM(1,8) ,NEXTST(256) ,THETA(64) ,KH(1,8) 
OPEN(UNIT=l , FILE= ' ERR ' , STATUS= ' OLD ' ) 

OPEN(UNIT=15,FILE='ERDATA' , STATUS= ' OLD ' ) 

WRITE (6,*) ' IF YOU WANT TO HAVE UNDETECTABLE ERROR' 

WRITE (6,*) 'ANALYSES ENTER "1" ,OTHERWISE ENTER "2"' 

READ (6,*) LD 

IF ( ORDER. EQ. 8) THEN 

KD=3 

ELSE 

KD=4 

ENDIF 

IF(LD.EQ.1)THEN 

WRITE (6,*)' WHAT KIND OF ANALYSES OUTPUT DO YOU WANT' 

WRITE (6,*)' (1) IN BINARY FORM' 



94 



13 

12 

333 

222 

444 

555 

45 

85 

55 

95 

10 



WRITE (6,*)' (2) IN DECIMAL FORM' 

READ (6,’^') IS 
WRITE (1,13) 

FORMAT(////,72{'#' ),//,7X, 'CORRECT RECEPTION COMBINATIONS IN PRESE 
*NCE OF ERROR' ,//,72( '#' )) 

WRITE (1,12) 

FORMAT (//,T5, 'INPUT' ,T17, 'PREVIOUS' ,T30, 'PRESENT' ,T43, 'PREVIOUS' , 
’^T55, 'PRESENT' ,T68, 'OUTPUT' ,/,T15, 'TRANSMISSION' ,T28, 'TRANSMISSION' 
’^',142, 'RECEPTION' ,T54, 'RECEPTION' ,/,79( ' -' )) 

ELSE 

GOTO 111 
END IF 
FL=2 

CALL XMTR (ORDER, FL) 

CALL XMSTA { THETA, TM, PS , ORDER) 

CALL :^-IDAT ( TM , ORDER , THETA , NEXTST , PS , T , TD ) 

DO 1 1=0, ORDER- 1 
DO 2 K=0, ORDER- 1 
A=TM(T(I) ,TD(K)) 

B=RCV1-IAT(TR(I) ,TR(A)) 

IF(IS.EQ.1)THEN 

OQU=A 

DO 333 JJ=KD,1,-1 

MOD(QQU,2) 



OU=K 

DO 222 LM=KD,1,-1 
KM(l,LM)=MOD(QU,2) 
QU=QU/2 
CONTINUE 
SQU=I 

DO 444 LN=KD,1,-1 
KN(l,LN)=MOD(SQU,2) 
SQU=SQU/2 
CONTINUE 



NT(1,JJ)= 

QQU=QQU/2 

CONTINUE 



SQU=B 

DO 555 LN=KD,1,-1 
KB(1 ,LN)=MOD(SQU,2) 

SQU=SOU/2 

continOe 

IF (ORDER. EQ.8)THEN 

WRITE (1, 45) (KM(1, L) ,L=1,KD), (KN(1,KK) ,KK=1 ,KD) , (NT(1 , IT) , IT=1 ,KD) 
*, (KN(1,KK),KK=1,KD) , (NT(1,IT) ,IT=1,KD) , (KB(l , IJ) , U=1 ,KD) 

FORMAT (T3, 3 (13) ,T15 ,3(13) ,T28,3(I3) ,T42,3(I3) ,T54.3(I3) ,T65,3(I3)) 
WRITEU5,85)(KN(1,KK),KK=1,KD),(NT(1,IT),IT=1,KD^ 
*,(KB(1,IJ).IJ=1,KD) 

FORMAT (9 (12)) 

ELSE 

WRITE(1,55)(KM(1,L),L=1,KD),(KN(1,KK),KK=1,KD),(NT(1,IT),IT=1,KD) 
*, (KN(1,KK) ,KK=1,KD) , (NT(1,IT) ,IT=1,KD) , (KB(1 , IJ) , U=1 .KD) 
F0RMAT(T1,4(I2),T15,4(I2),T28,4(I2) ,T42,4(I2) ,T54,4(I2) ,T65,4(I2)) 
WRITEU5,95)(KN(1,KK),KK=1,KD),(NT(1,IT),IT=1,KdJ 
*,(KB(1,IJ),IJ=1,KD) 

FORMAT (12 (12)) 

ENDIF 

ELSE 

WRITE(1,10) K,I,A,I,A,B 

FORMAT(T10,I3,T25,I3,T35,I3,T48,I3,T58,I3,T71,I3) 

ENDIF 

DO 3 L=0, ORDER- 1 
DO 4 J=0, ORDER- 1 
IF(L.EQ.I.AND. J.EQ.A)THEN 
GOTO 100 



ELSE 

Y=RCVMAT(TR(L) ,TR(J)) 
IF(Y.EQ.B)THEN 

if(is.eq.i)then 

SQU=L 

DO 666 LN=KD,1,-1 
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KH(1,LN)=M0D(SQU,2) 

S0U=SQU/2 
666 COfJTINUE 
SOU=J 

DO 777 LN=KD,1,-1 
KJ(l,LN)=MOD(SQU,2) 

SQU=SOU/2 
111 CONTINUE 
SOU=Y 

DO 888 LN=KD,1,-1 
KY{l,LN)=MOD(SQU,2) 

SQU=SOU/2 
888 CONTINUE 

IF(ORDER.EQ.8)THEN 

WRITE (1, 48) (KHU,KK) ,KK=1,KD) , (KJ(1,KK) ,KK=1,KD) 
(KY(1,IT) ,IT=1,KD) 

48 FORMAT (T42, 3(13) ,T54, 3(13) ,T65, 3(13) ) 

WRITE (15, 999 )(KH(1,KK) ,KK=1,KD) , (KJ(1,KK) ,KK=1,KD) 
Uy(1,IT) ,IT=1,KD) 

999 FORMAT (9(12)) 

ELSE 

WRITE(1,49)(KH(1,KK) ,KK=1,KD) , (KJ(1,KK) ,KK=1,KD) 
*,(KY(1,IT),IT=1.KD) ^ ^ ^ 

49 FORMAT (T42, 4(12) ,T54, 4(12) ,T65, 4(12) ) 

WRITE(15,949) (KH(1,KK) ,KK=1,KD) , (KJ(1,KK) ,KK=1,KD) 
*,(KY(1,IT),IT=1,KD) 

949 FORMAT (12(12)) 

END IF 
ELSE 

WRITE(1,20) L,J,Y 
20 FORMAT(T48,I3,T58,I3,T71,I3) 

ENDIF 

ENDIF 

100 ENDIF 
4 CONTINUE 
3 CONTINUE 

WRITE (1,46) 

46 FORMAT (79('-')) 

2 CONTINUE 
1 CONTINUE 
111 CLOSE(l) 

CLOSE (15) 

RETURN 

END 
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APPENDIX B 

COMPUTER GENERATED ENCODER AND DECODER OUTPUTS 
1. 8-DPSK TRANSMITTER 



TRANSMITTER'S ENCODER MATRIX TABLE FOR 8-DPSK 
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■> 



8-DPSK RECEIVER 



RECEIVER'S DECODER MATRIX TABLE FOR 8-DPSK 



GENERAL RECEIVER MATRIX 



01234567 

70123456 

67012345 

56701234 

45670123 

34567012 

23456701 

12345670 



RECEIVER MATRIX IS 
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J, 



16-DPSK TRANSMITTER 



TRANSMITTER'S ENCODER MATRIX TABLE FOR 8-DPSK 



GENERAL TRANSMITTER MATRIX 
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4, 



16-DPSK RECEIVER 



RECEIVER'S DECODER MATRIX TABLE FOR 16-DPSK 



GENERAL RECEIVER MATRIX 



012345678 9101112131415 
15 012345678 91011121314 
1415 012345678 910111213 
131415 012345678 9101112 
12131415 012345678 91011 
1112131415 012345678 910 
101112131415 0123456789 
9101112131415 012345678 
8 9101112131415 01234567 
7 8 9101112131415 0123456 
678 9101112131415 012345 
5678 9101112131415 01234 
45678 9101112131415 0123 
3 6 7 8 9101112131415 0 1 2 

2345678 9101112131415 0 1 
12345678 9101112131415 0 



RECEIVER MATRIX IS 



SYMBOL STATE (COLUMNS) / SYMBOL STATE (ROWS) 
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TRANSMITTER INFORMATION FOR THE GIVEN RECEIVER 



GENERAL TRANSMITTER MATRIX 



ABCDEFGHIJKLMNOP 

BCDEFGHIJKLMNOPA 

CDEFGHIJKLIWOPAB 

DEFGHIJKLMNOPABC 

EFGHIJKLMNOPABCD 

FGHIJKLMNOPABCDE 

GHIJKLMNOPABCDEF 

HIJKLMNOPABCDEFG 

IJKLMNOPABCDEFGH 

JKLMNOPABCDEFGHI 

KLMNOPABCDEFGHIJ 

LMNOPABCDEFGHIJK 

HNOPABCDEFGHIJKL 

NOPABCDEFGHIJKLM 

OPABCDEFGHIJKLMN 

PABCDEFGHIJKLMNO 



TRANSMITTER MATRIX 



SYMBOL STATE (COLUMNS) / DATA STATE (ROWS) 
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APPENDIX C 

COMPUTER GENERATED OUTPUTS FOR DBDR'S ERROR PATTERN 

ANALYSES 

The program outputs given below are a partial sample of the main outputs. 
Since undetectable error patterns are similar, in the next sections, 1/8 period of the 
pattern have been included. Actually 526 combinations of error patterns are detected 
by the program for the 8-DPSK case. For the 16-DPSK case there are 4096 
undetectable error pattern combinations. 

1. 8-DPSK RECEIVER 

CORRECT RECEPTION COMBINATIONS IN PRESENCE OF ERROR 



INPUT PREVIOUS PRESENT PREVIOUS PRESENT OUTPUT 

TRANSMISSION TRANSMISSION RECEPTION RECEPTION 
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I6-DPSK RECEIVER 



CORRECT RECEPTION COMBINATIONS IN PRECENSE OF ERROR 
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0 0 11 
110 1 
10 11 
10 0 1 
0 110 
10 0 0 
0 111 
110 0 
0 10 0 
10 10 
0 0 0 0 
0 0 10 



0 0 11 
1110 
0 0 0 1 
110 0 
10 11 
1111 
10 0 1 
110 1 
0 111 
10 10 



10 11 
10 11 
10 11 
10 11 
10 11 
10 11 
10 11 
10 11 



110 0 
110 0 
110 0 
110 0 
110 0 
110 0 
110 0 
110 0 
110 0 
110 0 
110 0 
110 0 
110 0 
110 0 
110 0 
110 0 



110 1 
110 1 
110 1 
110 1 
110 1 
110 1 
110 1 
110 1 
110 1 
110 1 
110 1 
110 1 
110 1 
110 1 
110 1 
110 1 



1110 

1110 

1110 

1110 

1110 

1110 

1110 

1110 

1110 

1110 

1110 

1110 

1110 

1110 

1110 

1110 



1111 

1111 

1111 

1111 

1111 

1111 

1111 

1111 

1111 

1111 
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10 10 
10 11 
.110 0 
•110 1 
1110 
1111 



0101 1111 
1000 1111 
0110 1111 
0010 1111 
0100 1111 
0000 1111 
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